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I. Wprowadzenie 


Praktyka dydaktyczna na Wydziale Energetyki i Paliw Akademii Górniczo-Hutniczej 
pokazała, że warto, aby studenci na początku studiów pogłębili zrozumienie drugiej zasady 
termodynamiki, nauczyli efektywnie posługiwać się takimi pojęciami, jak entropia, entalpia 
i entalpia swobodna. Niniejsza publikacja poszerza zakres materiału z przedmiotu fizyka II 
wykładanego na pierwszym roku studiów na uczelniach technicznych, daje wprowadzenie 
do termodynamiki maszyn przepływowych, takich jak turbina i kocioł, oraz termodynamiki 
chemicznej na przykładzie termodynamicznego opisu ogniwa paliwowego i syntezy czą- 
steczki amoniaku. 


I.1. Wymagania wstępne 


Przydatna jest znajomość podstaw termodynamiki, pierwszej zasady termodynamiki 
i własności gazu doskonałego w zakresie wymaganym na maturze z fizyki. Zagadnienia te 
są wyłożone w akademickich podręcznikach dla studentów pierwszego roku do przedmio- 
tu fizyka II, w tym podręczniku Fizyka dla szkół wyższych. Tom 2! (cz. 1: Termodynamika). 

Najważniejsze jest opanowanie materiału z trzech pierwszych rozdziałów, w których 
przedstawione są zagadnienia wstępne, własności gazu doskonałego i pierwsza zasada ter- 
modynamiki zapisana w konwencji AU = Q — W, gdzie W jest pracą wykonaną przez układ. 
Do właściwego zrozumienia prezentowanych zagadnień ważne jest również opanowanie 
chemii na poziomie podstawowym. 


I.2. Prosty silnik tłokowy podnoszący żwir 


Po opanowaniu podstaw można przystąpić do opisu prostego silnika tłokowego, w któ- 
rym gaz doskonały jest zamknięty w cylindrze z ruchomym, nieważkim tłokiem, poruszają- 
cym się bez oporów. Cylinder i tłok mają zerową masę, zerową pojemność cieplną, nie pod- 
legają odkształceniom, a gaz w cylindrze jest idealnie izolowany termicznie od otoczenia. 
Konstrukcja silnika pozwala zarówno podgrzewać, jak i chłodzić zawarty w nim gaz. Usta- 
wiając silnik pionowo, można go wykorzystać do podniesienia użytecznej masy na zadaną 
wysokość. Wybrano żwir ze względu na to, aby załadunek i rozładunek użytecznej masy 


: Fizyka dla szkół wyższych. Tom 2, Katalyst Education, OpenStax. Pobierz za darmo ze strony 


https://openstax.org/details/books/fizyka-dla-szkół-wyższych-tom-2-polska [dostęp: 15.11.2022] 


odbywał się kamyczek po kamyczku, by w każdej chwili panowały warunki zbliżone do 
stanu równowagi. 

W skrócie opisany silnik to idealny cylinder z ruchomym tłokiem, w którym zamknięto 
pewną ilość jednoatomowego gazu doskonałego i wiadomo, że w początkowym położeniu 
tłoka, w punkcie A pokazanym z lewej strony na rysunku I.1 parametry gazu są następujące: 


p,=p,=lbar=10Pa, V,=10dm*=0,01m, 7,=300K. 


W kolejnych charakterystycznych punktach pracy silnika, B, C, D, parametry gazu, 
położenie tłoka i jego obciążenie żwirem są takie, jak pokazano na rysunku 1.1. 


my = 0 kg mp = 100 kg mc = 100 kg Mp = 0 kg 
ha =fo =1m hp = ho hc = 2- hg hp =2-ho © 
A = 100 cm? A = 100 cm? A = 100 cm? A = 100 cm? 
V=% = 0,01m? V =V Ve =2+V Vp =2-V 4 
-p —105 
pa = Po = 10° Pa Ps = 2' Po Pc = 2: Po Pp = Po Ń 
Ta =T, =300K Ts =2:To és Tę = 4-Ty Tp =2:Ty Q 
AŻ 
bo o 2 g 
> © o O 
© v a Ee 
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oO = U = CcC 
w = ce co 
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Pat = Po = 10° Pa 


Rysunek I.1. Prosty silnik tłokowy podnoszący żwir 


Pomiędzy kolejnymi położeniami tłoka zachodzą procesy zaczynające się od grzania 
gazu i załadunku żwiru z zachowaniem stałej objętości (izochoryczny proces AB). W kolej- 
nych procesach kontynuowane jest ogrzewanie gazu w celu podniesienia obciążonego tłoka 
na zadaną wysokość (izobaryczny proces BC), rozładowanie żwiru na zadanej wysokości przy 
równoczesnym chłodzeniu gazu (izochoryczny proces CD) i dalsze chłodzenie gazu prowa- 
dzące do wepchnięcia tłoka przez ciśnienie atmosferyczne na początkową wysokość (izoba- 
ryczny proces DA). Zakłada się, że w każdym punkcie cyklu, czyli zarówno w punktach A, B, 
C, D, jak i w każdym punkcie pośrednim, spełniony jest warunek równowagi termodynamicz- 
nej, czyli w szczególności gaz w cylindrze zawsze ma dobrze określoną temperaturę, której 
wartość jest taka sama w całej objętości. Dalej zakłada się, że wszystkie procesy zachodzą 
w sposób rownowagowy’, czyli w każdej chwili wypadkowa siła działająca na tłok wynosi 


p 


W cytowanym podręczniku proces taki nazywany jest kwazistatycznym. 


zero. Taki statyczny, rownowagowy obraz stoi w sprzeczności z intuicyjnie postrzeganą dyna- 
miczną pracą silnika. Z całą pewnością trudno odrzucić takie intuicyjne wyobrażenie pracy 
silnika, ale jest to konieczne dla opanowania dalszego materiału. Można w tym celu wyobrazić 
sobie, że tłok w silniku porusza się bardzo powoli i że załadunek i rozładunek żwiru też prze- 
biegają bardzo wolno. Warto jednak dodać, że obraz powolnej pracy silnika będzie pomocny 
jedynie na obecnym, czyli początkowym etapie zgłębiania wiedzy o termodynamice. 

Zakłada się również, że przedstawiony silnik może pracować w obu kierunkach — jest 
maszyną odwracalną. Na początku przyjmijmy robocze, intuicyjne założenie, że odwracal- 
ność oznacza, że jeśli zostanie nakręcony film obrazujący cykl z podnoszeniem żwiru, to 
odtworzony od końca będzie opisywał silnik wykorzystany do opuszczenia żwiru. Jest to 
wyidealizowana sytuacja, gdyż zakłada zerowe opory ruchu tłoka oraz to, że warstwa usypa- 
nego na tłoku żwiru ma zerową wysokość. Niemniej wydaje się, że możliwe jest skonstruowa- 
nie rzeczywistego silnika, który będzie zbliżony do ideału, a film przedstawiający jego pracę 
puszczony od końca nie będzie raził ani śmieszył. Inaczej dzieje się w odniesieniu do proce- 
sów, które są w tak dużym stopniu nieodwracalne, że widać to gołym okiem. Na przykład, 
gdy samochód hamuje, to wiadomo, że rozgrzewają się jego tarcze hamulcowe. Film pusz- 
czony od tyłu pokaże rozpędzający się samochód i chłodzące się tarcze hamulcowe, co każdy 
uzna za coś śmiesznego i niemożliwego. Przykład ten dobrze ilustruje intuicyjne rozróżnienie 
procesów odwracalnych i nieodwracalnych, a w kolejnych rozdziałach niniejszej publikacji 
zostanie pokazane, w jaki sposób termodynamika pozwala na ilościowe ich rozróżnienie. 

Procesy zachodzące w opisanym idealnym, odwracalnym silniku można nakreślić na 
wykresie w przestrzeni ciśnienie-objętość, p-V (rys. 1.2). 


Vo 2Vo V 


AB — izochoryczne podgrzanie gazu do ciśnienia 2p, 
BC — izobaryczne rozprężenie do objętości 2V, 

CD - izochoryczne ochłodzenie do ciśnienia p, 

DA — izobaryczne sprężenie do objętości V, 


Rysunek I.2. Cykl w przestrzeni p-V prostego silnika ttokowego podnoszącego żwir 


Kolejne kroki w opisie działania silnika polegają na wykonaniu obliczeń dla tego cyklu, 
umieszczeniu wyników w tabelach oraz obliczeniu sprawności silnika (cyklu). 


1.2.1. Wykonanie obliczen dla prostego silnika tlokowego 
podnoszacego zwir 


Opisując silnik, podano ciśnienie, temperaturę i objętość gazu w punkcie A. Znając 
równanie stanu gazu doskonałego (równanie Clapeyrona) pV = nRT, można obliczyć liczbę 
moli gazu zgromadzonego w cylindrze. Dla wygody rachunkowej warto obliczyć iloczyn 
liczby moli i stałej gazowej: 

5 -2 
„ko PP OO N (1.1) 
Ti 3-107 3 

Podstawiajac dane liczbowe, nalezy pamietaé o odpowiednim przeliczeniu jednostek, 
gdyz w zadaniach czesto cisnienie jest wyrazane w barach (1 bar = 10° Pa), objetosci w decy- 
metrach sześciennych (1 dm? = 107 m*), a temperatura w stopniach Celsjusza, a nie w kel- 
winach, T,= Tot 273,15. 

Czynnikiem roboczym w silniku jest jednoatomowy gaz doskonały, którego molowe 
ciepło właściwe przy stałej objętości wynosi C, = (3/2)R, a przy stałym ciśnieniu C, = (5/2)R. 
Wartości liczbowe nC, i nC, odpowiednio wynoszą: 


5: 25 
nC, = nR- = 5 IK ka aa J/K (1.2) 


Ciśnienie w punktach B i C oraz w trakcie procesu BC jest sumą ciśnienia atmosferycznego 
i ciśnienia wynikającego z obciążenia tłoka o polu powierzchni A = 100 cm? = 10? m? żwirem 
o masie m = 100 kg. Przyjmując wartość przyśpieszenia ziemskiego g = 10 m/s”, dostaje się: 


+ 107-10 
Pp =P, ae aipa =2-10°Pa (1.3) 


Korzystając z podanych wartości i równania stanu gazu doskonałego (równania Cla- 
peyrona), można sprawdzić poprawność danych zawartych w tabeli I.1 parametrów gazu 
w punktach A, B, C, D. 


Tabela I.1. 
Parametry pracy silnika podnoszącego żwir w punktach A, B, C, D 


Punkt p [Pa] V [m°] T [K] 
A 1-105 1-10? 3:10? 
B 2-105 1-102 6: 107 
C 2-105 2-102 12:10? 
D 1-105 2-102 6: 107 


W izochorycznym procesie AB (V = const) gaz nie wykonał pracy, czyli pierwsza 
zasada termodynamiki o ogólnej postaci AU = Q — W przybiera postać: AU = Q,= nC AT. 
Po podstawieniu wartości liczbowych dostaje się: 


AU „=Q,=nC,AT,,= 5: 300 = 1500 J= 1,5 kJ (1.4) 


Czynnikiem roboczym silnika jest gaz doskonały, czyli zmiana energii wewnętrznej dla 
wszystkich procesów jest równa iloczynowi liczby moli, molowego ciepła właściwego przy 
stałej objętości i różnicy temperatur: 


AU,..=nC,AT,.= 5" 600 = 3000 J = 3,0 kJ (1.5) 
AU,,, = nC, AT. 5 : (600) = 3000 J = -3,0 kJ (1.6) 
AU, =nC,AT,, = 5 - (300) =-1500 J =—1,5 kJ (1.7) 


W izobarycznym procesie BC (p = const) ilość dostarczonego ciepła jest równa: 
Onc = NCAT, = (25/3) " 600 = 5000 J = 5,0 kJ (1.8) 


Teraz pracę wykonaną przez gaz można obliczyć, korzystając z pierwszej zasady ter- 
modynamiki: 


Wye = Qpe- AU ge = 5,0 — 3,0 = 2,0 kJ (1.9) 


lub z ogólnej zależności opisującej pracę wykonaną przez gaz zamknięty w cylindrze z rucho- 
mym tłokiem: 


Vo 
Wye = |p:dV = py AV,, = 2-10° -10° = 2000 J = 2,0 kJ (1.10) 


Vg 


albo wykorzystując równanie Clapeyrona, z którego wynika, że dla procesu izobarycznego: 


p: AV=nR- AT (1.11) 
czyli podstawiając: 


W,.=pz' AV c= nR: AT,.= (10/3) + 600 = 2000 J = 2,0 kJ (1.12) 


Izochoryczny proces CD przebiega analogicznie jak proces AB, ale zmiana energii 
wewnętrznej jest ujemna i ilość dostarczonego ciepła też jest ujemna. Podobnie izobaryczny 
proces DA przebiega jak proces BC, ale zmiana energii wewnętrznej jest ujemna, ponadto 
ujemną wartość ma ciepło dostarczone do gazu i gaz wykonuje ujemną pracę. Po wykonaniu 
niezbędnych prostych obliczeń dostaje się wartości zmian energii wewnętrznej AU, dostar- 
czonego ciepła O i wykonanej pracy W dla każdego z procesów AB, BC, CD i DA. Wyniki 
zapisano w tabeli I.2, aby ułatwić sprawdzenie ich poprawności i dalszą analizę. 


Tabela I.2. 


Zmiana energii wewnętrznej AU, dostarczonego ciepła Q i wykonanej pracy W 
w procesach AB, BC, CD i DA dla silnika podnoszącego żwir 


Symbol procesu AU [kJ] Q [kJ] W [kJ] 
AB 1,5 1,5 0,0 
BC 3,0 5,0 2,0 
CD —3,0 —3,0 0,0 
DA =S =25 =150 
Suma 0,0 1,0 1,0 
B FE =Q =W 

Interpretacja sumy suma musi być równa zeru na O * 

netto netto 
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Sumaryczna zmiana energii wewnętrznej po powrocie do punktu A musi być zerowa, 
gdyż podobnie jak wszystkie parametry stanu (temperatura, objętość, ciśnienie) również 
funkcja stanu, jaką jest energia wewnętrzna po wykonaniu pełnego cyklu, musi powrócić do 
wartości początkowych. Widać też, że praca netto musi być równa ciepłu netto, co wynika 
z pierwszej zasady termodynamiki, która przedstawia bilans energii. 


1.2.2. Obliczenie sprawności silnika podnoszącego żwir 


Z tabeli I.2 widać, że: 


O zostarczone ~ Lua = lQasl + Quel = 6,5 kJ (1.13) 

Om” Lac Lal "Ie = 5.5 kJ (1.14) 

Oeno = Das + Det Qen + Op, = 15 + 5,0 — 3,0 — 2,5 = 1,0 kJ (1.15) 

ann = Ogi = Oime z O wine = ,0kJ (1.16) 

veto 7 Wag t Wt Wen t Wpa = 9,0 + 2,0 + 0,0 — 1,0 =1,0 kJ (1.17) 
E Wae Oe ee „uj 


n= (1.18) 


Q, dostarczone LECIE O ae 


Warto dokładniej przeanalizować bilans pracy wykonanej przez gaz. W procesach AB 
i CD objętość gazu w cylindrze jest stała i gaz wykonuje zerową pracę. W procesie BC gaz 


wykonał dodatnią pracę: W,. = | p:dV = p, :AV,, =2 kJ. Graficznie w przestrzeni p-V 


powyższą całkę można pokazać jako pole ia pod odcinkiem BC, tak jak na 
rysunku I.3. 


X 
Vo WV 
Rysunek I.3. Proces BC podnoszenia tłoka ze żwirem w przestrzeni p-V 
Analogicznie w procesie DA gaz wykonał ujemną pracę, czyli zmniejszył swoją objętość: 


Va 
Wp = I p:dV =p,:NV,, = Pe" AV 5c =—1 kJ, co reprezentuje ujemne pole pod odcinkiem 
DA na rysunku I.4. 
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V WV 
Rysunek I.4. Proces DA opuszczania pustego tłoka w przestrzeni p-V 


Podsumowując, jak pokazano w równaniu (1.17), praca netto wykonana w jednym 
cyklu wynosi: W = 1 kJ i jest równa polu prostokąta ABCD w przestrzeni p-V, co ilu- 
struje rysunek 1.5. 


Vo WV 


Rysunek I.5. Praca netto cyklu ABCD w przestrzeni p-V jest rowna polu powierzchni prostokata 
ABCD, którego boki ilustrują przebieg kolejnych procesów cyklu 


Analizując dowolny cykl w przestrzeni p-V, można zauważyć, ze praca netto zawsze 
jest równa polu powierzchni figury zamkniętej i nakreślonej przez kolejne procesy cyklu. 
Jeśli kierunek kolejnych procesów cyklu obraca się w prawo, czyli zgodnie ze wskazówkami 
zegara, to praca netto jest dodatnia (w tym silniku, gdy podnosi on żwir). 

Na rysunku I.5 usiłowano również przedstawić przepływy ciepła dostarczonego Q,,, 
i odebranego Q, „p jednak niemożliwe jest zwłaszcza dokonanie wizualnego oszacowania 
ilości dostarczonego i odebranego ciepła, a tym samym wizualnego ocenienia sprawności 
tego cyklu, która — jak wynika z przeprowadzonych powyżej rachunków — jest mała, bo 
wynosi niewiele ponad 15%. 


II. Entropia 


II.1. S-— nowa wielkość charakteryzująca stan układu 


Bez wątpienia do oceny sprawności silników cieplnych przyda się wykres, na którym 
wprost będzie widać ilość ciepła dostarczonego do układu i odebranego z układu. W analizo- 
wanym przypadku układem jest gaz zamknięty w idealnym cylindrze. 

Aby taki rysunek przygotować, należy zauważyć, że ilość ciepła będzie zależała od tem- 
peratury układu, czyli w tym przypadku temperatury gazu w cylindrze. Dalsze rozważania 
prowadzą do wniosku, że należy wprowadzić nowy parametr opisujący własność gazu, który 
zostanie oznaczony literą S. Pole powierzchni pod funkcją T(S) od początkowej wartości 5, 
do końcowej, większej wartości S,, dla dowolnego odwracalnego i równowagowego procesu 
będzie równe ilości ciepła dostarczonego do układu, co pokazano na rysunku II.1. 


T 


Rysunek II.1. Ciepło jako pole pod krzywą w układzie współrzędnych 7-S 


Tym samym dostarczone do układu ciepło zostanie opisane całką: 


QO, = |r -dS (11.1) 
która pokazuje definicję wielkości S dla procesów równowagowych i odwracalnych: 
def 
ds E — (11.2) 


Wielkość S, nazywana entropia, jest funkcją stanu, gdyż charakteryzuje każdy rownowago- 
wy stan układu niezależnie od sposobu, w jaki został on osiągnięty. Identyczną własność mają 
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parametry stanu, czyli temperatura, ciśnienie i objętość, oraz poznana już energia wewnętrzna 
układu, która jest funkcją stanu. Różnica pomiędzy parametrem stanu i funkcją stanu dla przed- 
stawionych tu rozważań jest mało istotna i jest bardzo intuicyjna. W pewnym uproszczeniu 
parametrami stanu są wielkości bezpośrednio mierzone, takie jak temperatura, ciśnienie i obję- 
tość. Natomiast funkcjami stanu są wielkości obliczane na podstawie pomiarów, takie jak pozna- 
na wcześniej energia wewnętrzna i teraz wprowadzana entropia. Energia wewnętrzna gazu 
doskonałego zależy od temperatury 7, liczby moli gazu n oraz molowego ciepła właściwego C, 
i wyraża się prostym i dobrze już znanym wzorem: U = nC T. Co więcej, energia wewnętrzna 
ma też bardzo dobrą, intuicyjnie zrozumiałą interpretację nie tylko dla gazu doskonałego. Nato- 
miast zasadnicza trudność w zrozumieniu, czym jest entropia, wynika z dwóch faktów: 

1. Brakuje jakiejkolwiek naturalnej intuicji mówiącej, czym jest entropia. Pozostaje jedynie 
oprzeć się na użyteczności tego pojęcia, na obrazie w przestrzeni 7-5, na którym widać dla 
procesów równowagowych i odwracalnych ilość dostarczonego ciepła do układu jako pole 
pod krzywą opisującą przebieg odwracalnego, równowagowego procesu. 

2. Brakuje narzędzi do intuicyjnego szacowania entropii, takich jak termometr do pomiaru 
temperatury, która bardzo dużo mówi o energii wewnętrznej. Brakuje również uniwer- 
salnego aparatu matematycznego do obliczenia wartości entropii. 


Pokonanie pierwszej trudności, czyli intuicyjnego zrozumienia, czym jest entropia, war- 
to odłożyć na później — na okres po dobrym opanowaniu termodynamiki. Druga trudność 
w standardowych wykładach z termodynamiki zostaje pokonana przez wprowadzenie spój- 
nego, ale dość trudnego aparatu matematycznego opartego na analizie funkcji wielu zmien- 
nych i cząstkowych równaniach różniczkowych”. Okazuje się jednak, że dla procesów izo- 
termicznych (T = const) oraz dla gazu doskonałego o ustalonej liczbie moli można już teraz 
pokazać sposób obliczania zmian entropii. 


Dygresja: Należy rozważyć też drugą opcję, w której ciepło będzie równe polu 
powierzchnipodfunkcjąS(T)w granicach odpoczątkowej wartości temperatury 7, do 
końcowej, większej wartości 7,, czyli ciepło Q,, w równowagowym, odwracalnym 
procesie byłoby równe całce Q, = | S-dT. Przy takim założeniu widać jednak, 


że niemożliwe jest przekazanie ciepła w procesie izotermicznym. Tymczasem pro- 
cesy izotermiczne są podstawą termodynamiki, bez nich nie można wprowadzić 
cyklu Carnota, co tym samym wyklucza możliwość przyjęcia tej drugiej opcji. 
Warto w tym miejscu dodać, że proces izotermiczny jest jedynie wyidealizowaną 
koncepcją sprzeczną z obserwacją, że przepływ ciepła pomiędzy dwoma ciałami 
jest możliwy jedynie wtedy, gdy mają one różne temperatury. Aby pogodzić ze 
sobą te dwa wzajemnie wykluczające się obrazy, można przyjąć, że proces izo- 
termiczny to taki proces przepływu ciepła, w którym różnica temperatur jednak 
występuje, ale jest ona zaniedbywalnie mała. Widać jednak, że taki proces prze- 
wodzenia ciepła będzie nieskończenie powolny, co w jakimś sensie tłumaczy, dla- 
czego niemożliwe jest zbudowanie użytecznego silnika Carnota. Koniec dygresji. 


Warte polecenia są na przykład wykłady wideo i bardzo obszerne materiały pomocnicze bezpłat- 
nie udostępnione przez Massachusetts Institute of Technology (MIT): https://ocw.mit.edu/courses/ 
5-60-thermodynamics-kinetics-spring-2008/, a także MIT OpenCourseWare: https://ocw.mit.edu/. 
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II.2. Zmiana entropii dla odwracalnego, 
równowagowego procesu izotermicznego 


def 
Korzystając z definicji dS =6Q/T, można zauważyć, że dla odwracalnych, równo- 
wagowych procesów izotermicznych, czyli przebiegających przy stałej temperaturze 7, 
otrzymuje się: 
S, — S, = AS, = Lo (1.3) 
Na przykład” ilość ciepła potrzebna do stopienia kostki lodu o masie m w temperatu- 
rze 0° Celsjusza, Tœ wyniesie Q = m - C walić gdzie, © dzka jest ciepłem topnienia lodu. 
Widać zatem, że entropia układu (początkowo kostki lodu, a na końcu wody o temperatu- 
rze T c) zwiększy się o: 


m: C iopnienia 
AS (II.4) 


Tc 


Widać, że nawet w tak prostym zadaniu można obliczyć jedynie zmianę entropii układu, 
a nie bezwzględne jej wartości dla kostki lodu i wody. Okaże się, że zawsze obliczenie bez- 
względnej wartości entropii będzie albo bardzo trudne, albo wręcz niewykonalne. Na szczę- 
ście w zastosowaniach technicznych i chemicznych potrzebna będzie jedynie umiejętność 
obliczenia różnicy entropii pomiędzy dwoma stanami. 


II.3. Zmiana entropii gazu doskonałego 
dla dowolnego odwracalnego i rownowagowego procesu 


Pierwszą zasadę termodynamiki w ogólnym sformułowaniu AU = Q — W dla dwóch 
nieodległych stanów można zapisać jako: 


dU =80-8W (IL.5) 


Zastosowano tu konwencję zapisu, gdzie łacińska litera d oznacza różniczkę, grecka 
litera delta (ő) oznacza jedynie bardzo niewielką zmianę jakiejś wielkości. Na przykład jeśli 
w analizowanym procesie energia wewnętrzna układu zmieniła się o AU = 1000 J, to można 
przyjąć, że zmieniła się ona w milionie kroków, a każdy krok miał wartość 6U = 0,001 J. 
Natomiast symbol dU oznacza różniczkę, nieskończenie małą zmianę energii wewnętrznej, 
ale o dobrze określonych własnościach matematycznych, których omówienie wykracza poza 
ramy niniejszej publikacji. Należy jednak już teraz dodać, że dla takich wielkości jak ciepło 
i praca nie istnieją różniczki tych wielkości. Wprowadzona tu konwencja zapisu odzwiercie- 
dla fakt, że energia wewnętrzna jest funkcją stanu, natomiast ciepło i praca odnoszą się do 
procesów pomiędzy stanami i wielkości te zależą od metody dostarczenia ciepła i sposobu, 


4 


Patrz też przykład 4.4 w rozdziale 4.6: Entropia, w: Fizyka dla szkół wyższych. Tom 2, dz. cyt. 
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w jaki układ wykonuje pracę. Na szczęście dobre, matematyczne zrozumienie wprowadzonej 
konwencji nie będzie potrzebne w dalszej części niniejszej publikacji. Warto jednak już teraz 
stosować zaproponowaną konwencję, gdyż na dalszych etapach poznawania termodynamiki 
będzie to miało istotne znaczenie. 

W przypadku gazu zamkniętego w cylindrze z ruchomym tłokiem wiadomo, że dla 
stanów połączonych procesem odwracalnym i rownowagowym: dW = p : dV. Wprowadzając 
pojęcie entropii, dla rownowagowych i odwracalnych procesów otrzymano: Q = T : aS. 
Zatem dla takich procesów pierwszą zasadę termodynamiki opisuje równanie: 


dU=T:dS-p'dV (1.6) 


lub — gdy celem jest znalezienie entropii, czyli po podzieleniu równania (II.6) stronami 
przez T i po prostych przekształceniach — wzór: 


dS=—+4£.av (11.7) 


Jeżeli znane jest równanie stanu gazu, czyli zależność łącząca jego temperaturę, 
ciśnienie i objętość, to ze wzoru (II.7) można wyznaczyć funkcję opisującą zmiany entro- 
pii S. Najczęściej jest to zadanie trudne, gdyż już dla gazu van der Waalsa równanie sta- 
nu ma na tyle uwikłaną postać, że niezbędne jest wprowadzenie skomplikowanego apa- 
ratu matematycznego, co wykracza poza ramy niniejszej publikacji. Jeszcze większym 
wyzwaniem będzie znalezienie zmian entropii dla gazów rzeczywistych, w tym dla pary 
wodnej, co jest istotą termodynamiki technicznej. Jedynie dla gazu doskonałego równanie 
stanu, inaczej zwane równaniem Clapeyrona, ma na tyle prostą postać, pV = nRT, gdzie 
n oznacza liczbę moli gazu, R — stałą gazową, że możliwe jest już teraz wyznaczenie funk- 
cji opisującej zmiany entropii. Po pierwsze, z równania Clapeyrona widać, że zawsze speł- 
niona jest relacja: 


p nR 
—=— IL.8 
KICZ (11.8) 
Wiadomo też, że dla gazu doskonałego AU= nC, : AT, czyli dU=nC : dT. Podstawiając 
te dwie zależności do równania (II.7), dostaje się: 
aT dv 
+ 


dS =nC, : nR: (1.9) 
T V 


Jeśli liczba moli jest stała, to każdy człon tego równania zależy tylko od jednej funk- 
cji stanu lub jednego parametru stanu. Po lewej stronie równania jest entropia, po prawej 
w pierwszym członie jedynie temperatura, a w drugim członie tylko objętość. Zatem całko- 
wanie równania (II.9) od stanu początkowego p do stanu końcowego k można zapisać jako: 


SĘ Ty V, 
es = ae -— + ar — (11.10) 


gdzie w stanie początkowym temperatura, objętość i entropia wynoszą: T Vo Sp aw stanie 


. . p 
końcowym odpowiednio: 7,, V,, S,. 
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š 5 Xk n a yore 
Wiedzac, ze Í dx/x = In(x,/x „) oraz że niezerowa ilość gazu doskonałego zawsze ma 
Xp 


niezerową wartość temperatury, ciśnienia i objętości, można otrzymać: 


T, 
S,-S, ZAS = ct ol 7 (IL.11) 


p p 


Powyższe równanie to recepta na obliczenie różnicy entropii pomiędzy dowolny- 
mi dwoma stanami gazu doskonałego dla ustalonej liczby moli. Aby wykonać obliczenia, 
wystarczy znać parametry stanu, czyli liczbę moli, temperaturę, ciśnienie i objętość gazu. 
Korzystając z równania stanu pV = nRT, można zauważyć, że dla stałej liczby moli (np. dla 
gazu zamkniętego w cylindrze), dla dowolnych dwóch stanów zachodzi równość: 


Vie T 
Pate _ (1.12) 
PV, Ty 


Wykorzystując zależność (II.12), znając własności logarytmów i wiedząc, ze C=C,+R, 
można zapisać formułę (II. 11) również w postaci: 


V, 
S,-S, =AS,, =nC, -In| EE |+nC, n E (11.13) 
Pp Y, 
oraz: 
T, Pr 
S-S, AS, =nC, -In A —nR:ln| — (11.14) 
p P, 


W tym momencie warto przypomnieć adiabatyczny rownowagowy proces odwra- 
calny, w którym wymiana ciepła Q jest zerowa, czyli pierwsza zasada termodynamiki 
AU = Q — W ma postać: AU =—W. Korzystając z wprowadzonego powyżej pojęcia entropii, 
można pokazać, że dla takiego procesu zmiana entropii pomiędzy dowolnymi dwoma punk- 
tami będzie zerowa AS, = 0. 

„Powyższe wynika z definicji entropii dla odwracalnych, równowagowych procesów, 


dS = 6Q0/T, oraz z faktu, że gdy każdy elementarny przepływ ciepła 6Q na drodze od punk- 
tu początkowego p do punktu końcowego k jest zerowy, to w konsekwencji po zsumowaniu 
zmiana entropii pomiędzy stanami p i k też musi być zerowa AS „= 0. Korzystając z tej wła- 
sności odwracalnej adiabaty, można teraz wyprowadzić znane, szkolne równanie adiabaty: 
pV“ = const. Skoro AS, — 0, to z równania (II.13) wynika, ze 


y 
nC, -In| ŻE |+nC, -In| + |=0 (11.15) 
Pp V, 


def 
Po podzieleniu równania przez nC, i wprowadzeniu wyktadnika adiabaty x =C,,/C,: 
V, 
In| ŻE |+k-in| — |=0 (11.16) 
Pp y, 
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y K 
In| Pee -|=0 
PV, 
PE =P, 
czyli dla wszystkich punktów odwracalnej adiabaty zachodzi równość: 


pV"7 const (II.17) 


Można mieć nadzieję, że pokazana tu prostota wyprowadzenia równania adiabaty* choć 
częściowo rekompensuje trudy zapoznania się z entropią. 


II.4. Obraz cyklu w przestrzeni temperatura-entropia, T-S 


Można teraz obliczyć różnicę entropii pomiędzy dwoma dowolnymi stanami dla ide- 
alnego, odwracalnego silnika podnoszącego żwir na zadaną wysokość. W tym celu należy 
skorzystać z parametrów pracy rozważanego silnika podnoszącego żwir (tab. I.1) i obliczyć 
zmianę entropii względem dowolnie wybranego punktu, na przykład względem punktu A. 


AS, =nC, In| = |+nR-ln| = | = 3,5 J/K (11.18) 


V, 
AS jc = MC, -In| = |+ 7R- In| = |~ 9,3 JK 


AS, =nC, -In| 2 |+nR-ln ~5,8J/K 


Aby dokładnie narysować wykres dla tego silnika w przestrzeni 7-S, należy obliczyć 
zmiany entropii względem wybranego punktu A i wszystkich punktów cyklu, a nie jedy- 
nie punktów B, C i D. Korzystając choćby z arkusza kalkulacyjnego, można to wykonać 
z dowolnie dobrą dokładnością, czyli dla dowolnie dużej liczby punktów. Jednak już dyspo- 
nując wartościami dla punktów węzłowych, czyli znając AS» AS, AS',,, można naszkico- 
wać cykl w przestrzeni 7-S, łącząc punkty A, B, C, D odcinkami prostymi, co pokazano na 
rysunku II.2. 

Obraz cyklu w przestrzeni p-V pokazuje pracę wykonaną przez układ, a w przestrzeni T-S 
przepływy ciepła, co ilustruje rysunek II.2. W przypadku rozwazanego silnika podnoszące- 
go żwir widać, że jego sprawność jest bardzo mała. Praca netto, W = Deo = Quon ~ coa» 
równa jest polu powierzchni czworokąta ABCD w przestrzeni 7-S, czyli polu figury nakre- 
ślonej przez kolejne procesy cyklu. Pole to oznaczono kolorem zielonym i jest ono znacznie 
mniejsza od ilości dostarczonego ciepła Q,,,, wyróżnionego na czerwono. 


z Porównaj z rozdziałem 3.6: Proces adiabatyczny gazu doskonałego, w: Fizyka dla szkół wyższych. 


Tom 2, dz. cyt. 
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! Qnetto 
netto 


Sa SpSp Sc § 


Rysunek II.2. Praca jako pole pod krzywą w przestrzeni p-V 
i ciepło jako pole pod krzywą w przestrzeni T-S 


Warto też zobaczyć, jak wygląda cykl Carnota w przestrzeni T-S. W cyklu tym procesy 
są izotermami i adiabatami, czyli w przestrzeni T-S cykl Carnota ma kształt prostokąta, co 
pokazano na rysunku II.3. 


a) b) T 
Ty ot ~~ > m 
T cold PARKES 


A 'D 


Sa Sg ScsSp S 


Rysunek II.3. Cykl Carnota w przestrzeni p-V (a) i w przestrzeni 7-S (b) 
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W przestrzeni T-S, czyli na rysunku II.3, wprost widać, że wśród wszystkich cykli roz- 
piętych pomiędzy skrajnymi temperaturami 7,, i T. „ cykl Carnota ma największą spraw- 
ność, która wynosi: 


def = s 
n jal W eno = O esos Oddane = 1= T, Cold AS 4 = j= Too 
O esz PA T, Hot ` AS cg Tot 


Można mieć nadzieję, że pokazana tu prostota wyprowadzenia sprawności cyklu Carno- 
ta to kolejna zachęta do dokładniejszego poznania, czym jest entropia. 


(11.19) 


III. Idealny silnik cieplny, czyli jaki? 


W cytowanym już podręczniku” pod koniec rozdziału 4.4 Sformułowania drugiej zasa- 
dy termodynamiki pada stwierdzenie: 


Podsumowując, dochodzimy do wniosku, że wszystkie silniki z obiegiem odwra- 
calnym, działające między tymi samymi rezerwuarami, mają te same sprawności. 


Tymczasem, jak pokazano w poprzednich rozdziałach, sprawność silnika do podnosze- 
nia żwiru zasilanego przez dwa rezerwuary ciepła o skrajnych temperaturach T, = 1200 K 
oraz Toy = 300 K niewiele przekracza 15%. Natomiast gdyby pomiędzy tymi rezerwuara- 
mi ciepła pracował silnik Carnota, to jego sprawność wynosiłaby neo = 1— Teoy/T, 


= | — 300/1200 = 75%. Oczywiście oba rozważane silniki (cykle) są odwracalne, wice wade 
sprzeczność pomiędzy cytowanym z podręcznika zdaniem oraz przedstawionymi powyżej 
obliczeniami. Jest to jedynie pozorny paradoks wynikający z różnego określenia, co jest 
układem. W cytowanym powyżej rozdziale podręcznika układem (silnikiem) jest silnik wła- 
ściwy wraz z rezerwuarami ciepła zimnego i gorącego. Natomiast w rozważanym powyżej 
zadaniu z silnikiem podnoszącym żwir układem był sam silnik. Zagadnienie to z jednej stro- 
ny jest na tyle ważne dla zrozumienia drugiej zasady termodynamiki, że warto poświęcić mu 
więcej uwagi. Z drugiej strony jest na tyle trudne i czasochłonne, że jest jedynie materiałem 
dodatkowym dla osób pragnących poznać wyjaśnienie zagadki i bez utraty ciągłości prezen- 
towanego w publikacji materiału można teraz bezpośrednio przejść do rozdziału IV. 


III.1. Mikrosilnik Carnota 
pomiędzy dwoma zbiornikami ciepła 


Mikrosilnik Carnota o pomijalnie małej masie czynnika roboczego zostanie ustawio- 
ny pomiędzy dwoma zbiornikami wody, co pokazano na rysunku III.1. Zimny zbiornik ma 
masę m,i początkową temperaturę 7,,, a gorący ma masę m, i temperature 7,,. Oba zbior- 
niki i silnik są idealnie izolowane termicznie od otoczenia, są układem zamkniętym (mają 
stałą masę), ale układ ten wyprowadza na zewnątrz użyteczną pracę. Silnik jest pomijalnie 
mały i w każdym pojedynczym mikrocyklu temperatury zbiorników Hot i Cold są stałe, ale 
po każdym cyklu T, zmniejsza się o niewielką wartość ôT, a 7,. zwiększa o 67... Po pewnym 
czasie obydwa zbiorniki osiągną temperaturę równowagi 7, i silnik się zatrzyma. 


€ Fizyka dla szkół wyższych. Tom 2, dz. cyt. 
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Rysunek II.1. Mikrosilnik Carnota pracujący pomiędzy dwoma zbiornikami ciepła 


Zewnętrzny obserwator nie wie, w jaki sposób praca jest generowana, i jedynie czerpie 
korzyści z wykonanej przez silnik i dwa zbiorniki ciepła pracy mechanicznej. Widzi też, że 
układ ma stałą masę oraz że wraz w upływem czasu temperatury zbiorników zasilających 
silnik umieszczony pomiędzy zbiornikami się wyrównują. 


III.1.1. Temperatura równowagi 


Założono, że silnik umieszczony pomiędzy zbiornikami pracuje w cyklu Carnota, czyli 
jego sprawność w każdym mikrocyklu można wyrazić, podając zarówno przepływy ciepła, 
jak i temperatury wody w zbiornikach: 


= 50, — 00, = T =o 
50, Ty 


gdzie Q, jest ilością ciepła dostarczonego do silnika, a 6Q, odebranego w jednym 
mikrocyklu pracy, gdy temperatury wody w zbiornikach dla tego mikrocyklu wynosiły 
odpowiednio 7, i To. Wzór (II.1) sprowadza się do równości 60./60,, = T/T, i należy 
pamiętać, że w przyjętej konwencji obie wielkości ö0,, i 6Q, są dodatnie. W cyklu Car- 
nota ciepło jest pobierane i oddawane w procesach izotermicznych. W tym przykładzie po 
każdym cyklu temperatury zbiorników zmieniają się o wartości odpowiednio wynoszące 
67,1 67,,. Układając odpowiednie równanie różnicowe, należy trzymać się ustalonej kon- 
wencji znaków. Wydaje się, że najwygodniej myśleć, że zarówno 67,, jak i 67, są wiel- 
kościami dodatnimi. Proces grzania i chłodzenia wody w zbiornikach opisują równania 
kalorymetryczne: 


(IIL.1) 


50, m mc ÒT, 
Q- = mc ÒT, 


(IIL.2) 


gdzie symbolami m, i m, oznaczono odpowiednio masy wody w gorącym i zimnym zbiorni- 
ku, ac, jest ciepłem właściwym wody, które z założenia jest wielkością stałą. 
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Po podstawieniu do wzoru 60./5Q,,= T,./T, dostaje się równanie różnicowe: 


my: OT, = T, 


c 
m IL. 
m,°8T, Ty „A 


Chcąc wykorzystać rachunek różniczkowy, należy przekształcić równanie różnicowe 
w równanie różniczkowe. Operacja ta sprowadza się do ustalenia znaku przy transpono- 
waniu dodatnich, małych wielkości òx na nieskończenie małe różniczki dx występujące 
w rachunku różniczkowo-całkowym, gdzie zawsze symbol dx oznacza, że od wartości póź- 
niejszej odejmuje się wartość wcześniejszą. Układając równanie różnicowe, przyjęto, że ôT, 
i ôT, są wielkościami dodatnimi, ale wiadomo, że temperatura zbiornika zimnego rośnie, 
czyli dodatnia zmiana temperatury ôT „ będzie oznaczała również dodatnią zmianę dT: 


81. +dT, (IIL.4) 


Natomiast temperatura zbiornika gorącego maleje, czyli d7,, jest wielkością ujemną. 
Aby òT, było dodatnie, należy przyjąć: 


òT, > —dT, I.5) 
Po tym podstawieniu równanie różnicowe przechodzi w równanie różniczkowe: 
m, dT, T, 
— 2 (111.6) 
m, aT, To 


Niektóre równania różniczkowe, na przykład (II.6), można rozwiązać, rozdzielając 
zmienne, czyli zapisując je w taki sposób, aby każdy jego człon zależał tylko od jednej 
zmiennej. Po prostych przekształceniach dostaje się równanie, w którym po lewej stronie jest 
jedynie zmieniająca się temperatura 7,., a po prawej temperatura T; 


m, aT. z _ aT, 
H To Ly 


(111.7) 


m 

Teraz obie strony równania można niezależnie zsumować, czyli scałkować w granicach 
od stanu początkowego do końcowego. Po lewej stronie stan początkowy opisany jest przez 
początkową temperaturę zbiornika Cold wynoszącą T a po prawej stronie równania zbiorni- 
ka Hot o temperaturze początkowej wynoszącej 7. Stan końcowy po obu stronach równania 
określa szukana temperatura równowagi wynosząca 7,. Równanie całkowe ma zatem postać: 


T, T, 
m. ¢ dT, t dT, 
L. f=—|—# (I.8) 
My Teo Te Tyo H 
na = 
Wiedząc, że Í “= in( dostaje sie: 
M 6 x, 
T, T, 
mc af R Jn R | (IIL.9) 
My Teo Tao 
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Z wykorzystaniem własności logarytmow można to zapisać jako: 


mę My 


mc+my me +my III. 10 
Tr = Too Tao ( ) 


Wprowadzenie różnych mas wody zimnej i gorącej służyło lepszemu rozróżnieniu tych 
dwóch zbiorników, ale w dalszej analizie prowadzi jedynie do komplikacji rachunkowych. 
Dlatego dalej rozważania będą prowadzone dla równych mas, czyli dla m „= m „= m. Widać, 
że w tym szczególnym przypadku dostaje się: 


Tr = 4TyoTco (III.11) 


Warto zauważyć, że temperatura równowagi 7, jest równa średniej geometrycznej 
początkowych temperatur zbiornika zimnego i gorącego. 


II.1.2. Analiza entropii 


Analiza entropii zostanie przeprowadzona dla całego systemu, ale najpierw dla każdego 
elementu oddzielnie, czyli dla silnika oraz dla zbiornika zimnego i gorącego. Entropia jest 
wielkością ekstensywną i całkowita jej zmiana to suma zmian entropii poszczególnych ele- 
mentów systemu: 


AS 


układu 


= AS 


silnika 


+ AS 


zbiornika Hot 


+ AS 


zbiornika Cold 


(IIl.12) 


Zmiana entropii silnika będzie zerowa, gdyż jest to mikrosilnik o zaniedbywalnie małej 
wielkości. 

Teraz uwaga zostanie skoncentrowana na zbiornikach, pomiędzy którymi następuje 
przepływ ciepła. Na początku zostanie przedstawiona analiza zmiany entropii dla sponta- 
nicznego przewodzenia ciepła. Innymi słowy w tym przypadku nie ma między zbiornika- 
mi żadnego silnika, a przepływ ciepła nastąpi drogą przewodzenia. Jest to spontaniczny 
proces wyrównywania temperatur, czyli proces nieodwracalny. Z prostych rozważań kalo- 
rymetrycznych dostaje się, że dla równych mas wody w zbiornikach temperatura rów- 
nowagi jest równa średniej arytmetycznej początkowych temperatur zbiornika gorącego 
i zimnego: 


T _ Tao +Tco 


R przewodzenie ~ 2 


(111.13) 


Dla każdego ze zbiorników oddzielnie proces pobierania i odpowiednio odbierania cie- 
pła jest opisany zależnością kalorymetryczną. Jeśli podzielić ten proces na wiele mikroproce- 
sów, to w przypadku zbiornika gorącego każdy z nich prowadzi do przepływu ciepła: 


8Q,„= mc ôT, (111.14) 


gdzie c, jest ciepłem właściwym wody. Co więcej, przepływ ciepła w każdym takim 
mikroprocesie jest odwracalny, gdyż bez żadnej wątpliwości niewielka różnica temperatur, 
ale o przeciwnym znaku spowoduje identyczny przepływ ciepła, ale w przeciwną stronę. 
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Czyli zmianę entropii zbiorników z wodą można obliczyć ze znanej zależności dla procesów 
odwracalnych: 


mc dT, mc, dT, 
dS „piornika Hot T a , dS „pornika Cold = = (IT. 15) 
H È 


Widać, że znak zmian entropii jest taki sam jak znak zmian temperatury, a całkowita 
zmiana entropii układu składającego się z dwóch zbiorników wyniesie: 


AS kladu przewodzenie = AS ika Hot F AS biornika Cold (III. 16) 

TR przewodzenie TĄ przewodzenie d T 
= mc — 4 mc Ś (01.17) 

układu przewodzenie w w T 

Tuo H Teo [63 

pe 
R R przewodzenia 1.18 
AS [adu przewodzenie mc, In T ( E ) 
HO CO 


gdzie: Tr przewodzenie = m + Too)/2, co daje: 


(III.19) 


(Tao + Too) | 
AT 01 co 


AS adu przewodzenie = mc, | 


Widać, że zmiana entropii układu dwóch zbiorników, gdy doszło do wyrównania 
wartości temperatur zbiorników Hot i Cold na drodze przewodzenia, będzie dodatnia, 
ukladu przewodzenie 7 0» CO wynika z faktu, że średnia arytmetyczna dwóch różnych dodatnich 
liczb jest zawsze większa od ich średniej geometrycznej (należy pamiętać, że aby silnik 
ruszył, to początkowe temperatury zbiorników Hot i Cold muszą być różne i musi zostać 
spełniony warunek 7,, > To). 

Oba zbiorniki były termicznie odizolowane od otoczenia, więc przepływ ciepła z oto- 
czenia do tych zbiorników był zerowy. Co więcej, jeśli proces wyrównywania temperatur 
zostanie podzielony na mikroprocesy, to dla każdego z nich 00 , , ua = 0. Zatem w przy- 
padku układu dwóch zbiorników dla każdego mikroprocesu prowadzącego do wyrównania 
temperatur przez przewodzenie ciepła zachodzi nierówność: 


50% układu przewodzenia (II 20) 


układu przewodzenie T 


Analiza zmiany entropii układu dwóch zbiorników potwierdziła, że spontaniczny pro- 
ces wyrównania temperatur pomiędzy zbiornikami przez przewodzenie ciepła jest procesem 
nieodwracalnym, gdyż entropia układu wzrosła, pomimo że do układu dostarczono zerową 
ilość ciepła w każdym z mikroprocesów składających się na proces wyrównania temperatur, 
czyli gdy dla każdego mikroprocesu Q , am 7 9 


Powyższa analiza pokazała również, że aby proces wyrównania temperatur pomiędzy 
zbiornikami o takiej samej masie był odwracalny, czyli o zerowej zmianie entropii ukła- 
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du dwóch zbiorników, to temperatura równowagi musi być równa średniej geometrycznej 
początkowych temperatur zbiorników: 


Hof co (III.21) 


W pierwszej części rozwiązywanego zadania, w równaniu (III.11) pokazano, że jeśli do 
wyrównania temperatur posłuży mikrosilnik Carnota zainstalowany pomiędzy zbiornikami, 
to zostanie osiągnięta właśnie taka temperatura równowagi. Czyli proces przeprowadzający 
układ składający się z dwóch zbiorników i mikrosilnika Carnota od stanu początkowego, do 
stanu temperatury równowagi jest procesem odwracalnym. 

Pozostaje zadać pytanie: co się stanie, jeśli pomiędzy zbiornikami będzie pracował inny 
odwracalny silnik, na przykład silnik Stirlinga? Wówczas na zewnątrz zostanie wyprowa- 
dzona mniejsza ilość pracy, bo każdy silnik inny niż silnik Carnota ma mniejszą sprawność. 
Oznacza to, że temperatura równowagi będzie wyższa niż ta osiągnięta za pomocą silnika 
Carnota, czyli dochodzenie do temperatury równowagi układu „silnik plus dwa zbiorniki”, 
będzie mieć charakter nieodwracalny. Innymi słowy powyższa analiza pokazuje, że silnik 
Carnota jest jedynym idealnym silnikiem wśród wszystkich odwracalnych silników. 


II.2. Silnik Carnota — jedyny idealny silnik 


Przedstawiona analiza umożliwia dobrą interpretację wielu fragmentów z podręczni- 
ków akademickich, w tym z cytowanej już Fizyki dla szkół wyższych. Tom 2. W największym 
stopniu dotyczy to rozdziału 4.4: Sformułowania drugiej zasady termodynamiki. Padają tam 
też takie zdania, jak: 


Druga właściwość jest następująca: wszystkie silniki o obiegu odwracalnym, dzia- 
łające między dwoma tymi samymi rezerwuarami, mają taką samą sprawność. [...| 


Teraz dość łatwo jest udowodnić, że sprawności wszystkich silników z obiegiem 
odwracalnym, działających między tymi samymi rezerwuarami, są równe. 


Należy pamiętać, że są one poprawne przy założeniu, że silnikiem jest układ składający 
się z silnika właściwego i zbiorników ciepła. Wówczas, tak jak pokazano w rozdziale III. 1, 
dowodzi się, że istnieje tylko jeden idealny, odwracalny silnik i jest nim silnik Carnota. 

Porównując zmiany entropii pomiędzy dwoma dowolnie bliskimi stanami układu, moż- 
na ilościowo rozróżnić procesy odwracalne, dla których dS = 6Q/T, od nieodwracalnych, gdy 
dS > 6Q/T, co stanowi istotę drugiej zasady termodynamiki. Jeśli badana jest odwracalność 
procesów łączących odległe stany A i B, to drogę łączącą te stany należy podzielić na nie- 
skończenie małe odcinki. Jeśli dla choć jednego z nich zajdzie nierówność dS > 60/T, to cały 
proces jest nieodwracalny. Jeśli jednak wiadomo, że w trakcie procesu od stanu A do B nie 
doszło do wymiany ciepła, to badając odwracalność procesu AB, wystarczy sprawdzić, czy 
AS ,, = 9. Stosując tę metodę, należy jednak się upewnić, czy nie doszło do żadnej wymiany 
ciepła. Sam zerowy bilans przepływu ciepła na drodze od stanu A do stanu B to za mało, gdyż 
na przykład w początkowej fazie ciepło mogło być pobierane, a następnie oddawane. 
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IV. Idealne, odwracalne maszyny przeplywowe 


W standardowym programie przedmiotu fizyka II maszyny przepływowe nie są oma- 
wiane. Jednak specyfika Wydziału Energetyki i Paliw AGH wymaga choć wstępnego przed- 
stawienia termodynamiki idealnej turbiny, sprężarki, kotła i chłodnicy. 


IV.1. Idealna turbina i sprężarka 


Odizolowane termicznie od otoczenia idealna turbina i sprężarka są odwracalnymi, 
przepływowymi maszynami cieplnymi. Idealna turbina puszczona w odwrotnym kierun- 
ku jest idealną sprężarką o tych samych parametrach pracy. Przyjmując konwencję zgodną 
z pierwszą zasadą termodynamiki, turbina wykonuje dodatnią użyteczną pracę na zewnątrz 
układu, a sprężarka pracę ujemną, gdyż dodatnią pracę wykona siła zewnętrzna, na przykład 
silnik elektryczny ją napędzający. 

Rozważona zostanie idealna turbina zasilana gazem, przez którą przepływa ń moli gazu 
na sekundę i która pracuje ze stałą mocą oznaczaną jako P „„ co ilustruje rysunek IV.1. Każ- 
dy mol gazu, wpadając do turbiny, ma ciśnienie Ping? objętość molową (czyli objętość jednego 
mola gazu) v,,„» oraz temperaturę T, „a opuszczając turbinę ma parametry p „„ Vao Tue Należy 
podkreślić, że zostanie pominięta energia kinetyczna paczek gazu. Innymi słowy na począt- 
kowym etapie studiowania termodynamiki maszyn przepływowych zakłada się zerową pręd- 
kość unoszenia, z jaką gaz przepływa przez turbinę. 


P, turb [w] 


n [mol/s] n [mol/s] 


Pinp Vinp Pout Vout 
T; 7 


LI 
l 
1 
l 
U 


1 
I 
1 
1 
1 
1 


Zerowa wymiana u 
ciepła pomiędzy 
turbiną i otoczeniem 


Rysunek IV.1. Idealna turbina 


28 


Zatozenie to jest dobrze uzasadnione, gdyz wplyw energii kinetycznej kazdej pacz- 
ki z gazem jest niewielki, co można wykazać. Z drugiej strony założenie to jest trudne do 
zaakceptowania, gdyż przez turbinę przepływa ń moli gazu na sekundę. Oczywiście bar- 
dziej dokładne obliczenia będą wymagały uwzględnienia pominiętej tu energii kinetycznej 
przepływających paczek z gazem, ale to wymaga poznania podstaw termodynamiki maszyn 
przepływowych, co jest celem niniejszej publikacji. 

Ze schematu pokazanego na rysunku IV.1 widać, że idealna turbina pracuje w sta- 
bilnych warunkach — ani nie grzeje się, ani nie stygnie. Oznacza to, że energia zgroma- 
dzona w turbinie E „, (w jej korpusie, wirniku, gazie, który jest w środku, suma energii 
mechanicznej i ciepła) jest stała, czyli pochodna po czasie energii zgromadzonej w tur- 
binie jest zerowa: 


d 
— (E,„) = 0 IV.1 
aS sz) ( ) 


Do turbiny wpada moli gazu na sekundę, a każdy mol wnosi u, energii wewnętrznej. 
Należy jednak pamiętać, że gaz musi zostać wepchnięty do turbiny, co zwiększy energię 
w niej zgromadzoną o pracę objętościową, w, = Poin Analogicznie należy uwzględnić 
pracę wyrzucenia gazu z turbiny. 

Przyjmując konwencję, że turbina wykonuje dodatnią pracę na zewnątrz (np. napędza 
generator) ze stałą dodatnią mocą P, „, równanie (IV.1) bilansujące pochodną po czasie ener- 
gię zgromadzoną w turbinie można zapisać w postaci: 


inp 


d ; 
0 = dt (Enr) =A: (CH = Pimp Vinp ) m (U + Pow Vout » z Paro (IV.2) 


Widać, że warto wprowadzić nową wielkość, entalpię molową, czyli entalpię jednego 
mola gazu: 


def 


h=u+ pv (IV.3) 


def 
Analogicznie dla dowolnej ilości moli gazu definiowana jest entalpia: H =U + pV. 
Zastosowano konwencję oznaczeń, w której małymi literami oznaczono wielkości eksten- 
sywne dla jednego mola gazu i zachodzą równości: V = n: v, U=n:u, H- n : h oraz dla 
entropii S = n : s, gdzie n oznacza liczbę moli. Korzystając z pojęcia entalpii molowej i prze- 
kształcając wzór (IV.2), moc turbiny możemy zapisać: 
P =ń:(h 


turb 


h) (IV.4) 


inp = 

Aby zachować zgodność znaków z rachunkiem różniczkowym, gdzie zawsze od wiel- 
kości późniejszej odejmowana jest wcześniejsza, zmianę entalpii Ah zdefiniowano jako róż- 
nicę entalpii gazu wychodzącego i wchodzącego do turbiny: 


Ah=h —h (IV.5) 


out inp 


Daje to moc turbiny dla dowolnego gazu rzeczywistego: 


=- Ah (IV.6) 
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Skoro idealna sprezarka i idealna turbina sa wzgledem siebie odwracalnymi maszyna- 
mi, to moc sprężarki P „„ wyraża się tą samą zależnością: 


=P. 


spre turb 


-ñ Ah (IV.7) 


Innymi słowy idealna turbina od idealnej sprężarki różni się jedynie tym, ze moc turbiny 
jest liczbą dodatnią, gdyż temperatura gazu po przejściu przez turbinę zmniejsza się i AA jest 
ujemne, a moc sprężarki jest liczbą ujemną. 

Dla gazu doskonałego uzyskany wzór (IV.7) opisujący moc turbiny można dalej prze- 
kształcić, gdyż wówczas Au = C AT, gdzie AT= TT, — Lap Stosując równanie Clapeyrona dla 
jednego mola gazu, pv = RT, dostaje się: 


Ah = A(u + pv) = Au + A(pv) = CAT+ RAT = (C, + R) AT (IV.8) 


gdzie skorzystano z zależności A(pv) = RAT wynikającej z równania Clapeyrona. Dalej, 
wykorzystując równość C,=C,+R, dostaje się: 


Ah = CAT (IV.9) 
i moc idealnej turbiny zasilanej gazem doskonałym: 


P =ñ: CAT (IV.10) 


tur! 


W zastosowaniach technicznych przyjęto podawać przepływ masowy gazu m w kilo- 
gramach na sekundę. Ponieważ masa cząsteczkowa gazu M, najczęściej wyrażana jest w gra- 
mach, przepływ molowy i masowy łączy relacja: 


1000-75, 
5 _ g/s] 
Mi mois] e M (IV. 1 1) 


A [s] 


Idealna turbina jest maszyną odwracalną i izolowaną termicznie od otoczenia. Można 
zatem wnioskować, że stan wejściowy gazu Boa R musi leżeć na tej samej odwracalnej 
adiabacie co stan wyjściowy gazu opisany parametrami p „p V,,,, T „e Jeśli turbina napędzana 
jest gazem doskonałym, to parametry te muszą być związane znaną relacją By. T Pou Y oue 
gdzie k = C/C „jest wykładnikiem adiabaty. Jednak dla maszyn przepływowych objętość 
molowa v jest niewygodnym w użyciu parametrem i korzystając z funkcji stanu dla jednego 


mola, pv = RT, powyższy związek można wyrazić za pomocą ciśnienia i temperatury: 


x 1x 


P Typ = Pout To (IV. 12) 


Wzór (IV.12) umożliwia na przykład obliczenie b gdy znane są ciśnienia Die Pou 
i temperatura T, Ten bardzo użyteczny wzór trudno zapamiętać i warto poznać metodę jego 
wyprowadzenia bazującą na bilansie zmian entropii. Ogólnie metoda rozwiązywania zadań 
oparta na bilansie entropii jest w wielu przypadkach bardzo wydajna i choćby dlatego warto 
ją poznać, zaakceptować i używać. Podobnie jak w przypadku energii zgromadzonej w tur- 
binie, również entropia zgromadzona w turbinie, S, „, musi być stała, czyli jej pochodna po 
czasie musi być zerowa. 
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d 
— (S) = 0 IV.13 
ue) (IV. 13) 


Do turbiny wpada 7 moli gazu na sekundę: każdy mol wnosi > entropii, a wylatują- 
cy mol gazu wynosi s „ entropii. Z założenia, że turbina jest maszyną odwracalną, wynika 
jedynie, że jest dopuszczalna wymiana ciepła pomiędzy turbiną i otoczeniem, ale pod warun- 
kiem, że dla każdej porcji wymiany ciepła 0q, ran, pomiędzy turbiną i otoczeniem spełniona 
jest relacja charakteryzująca procesy odwracalne: 


O, 
wymiany 
W a (Tv. 14) 


wymiany 


Jednak dla idealnej turbiny założono, że nie ma wymiany ciepła pomiędzy turbiną i oto- 
czeniem i że w szczególności idealna turbina nie wymaga chłodzenia. Czyli całkowity brak 
wymiany ciepła z otoczeniem dla idealnej turbiny oznacza, że każda elementarna wymiana 
ciepła z otoczeniem jest zerowa, dq, „= 0, czyli zawsze ds, , = 0, a zatem całkowita 

. = 7 , Wymiany Boz: RY wymiany Š > 
wymiana entropii z otoczeniem też jest zerowa i jej pochodna po czasie też jest zerowa, 

wizę © Zatem ogólne równanie bilansujące pochodną po czasie entropii zgromadzonej 
w turbinie ma postać: 


d A 
0 = dt (Sr) = A . (Srp Sout ) E Bomian (IV. 15) 
ale jesli ee = 0, czyli dla idealnej, adiabatycznej turbiny sprowadza się do równania: 
out m S inp = As = 0 (IV. 16) 


Warunek ten ma zastosowanie w odniesieniu do kazdego gazu napedzajacego idealna, 
odwracalna, adiabatyczną turbinę. Jednak jedynie w przypadku gazu doskonałego różnicę 
entropii pomiędzy stanem końcowym i początkowym każdego mola gazu można wyrazić 
znaną relacją (II. 14), która dla jednego mola gazu, czyli gdy n = 1, ma postać: 


T, 
Ś _ Sinp = As = C, u out | R ol B= ) (IV. 17) 


inp P inp 


Jak łatwo sprawdzić, relacja (IV.17) jest tożsama z wcześniej otrzymanym wzorem 
l-K 1-« 


z równania adiabaty dla gazu doskonałego (IV. 12): BRL, =pxT, 


out * 


IV.2. Idealny kocioł i idealna chłodnica 


Podobnie jak turbina i sprężarka zostanie opisany idealny kocioł i idealna chłodnica, 
czyli druga para maszyn przepływowych wzajemnie odwracalnych. 

Rozważony zostanie idealny kocioł, przez który przepływa ri moli gazu na sekundę. Jest 
on podgrzewany ze stałą mocą O co ilustruje rysunek IV.2. Każdy mol gazu, wpadając 


grzania” 
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do kotła, ma ciśnienie Prp objętość molową V, Oraz temperaturę T, a opuszczając kocioł, 
ma parametry: p, = Pmp Your Tour Należy podkreślić, że podobnie jak dla turbiny pomijana 


jest energia kinetyczna paczek gazu. 


n [mol/s] n [mol/s] 
Pinp Vinp Pout Vout 
Tinp 1 \ out 

—> | => 


tttt Q grzania [W] 


Rysunek IV.2. Idealny kocioł 


. 
La 


Ze schematu widać, że podobnie jak w przypadku idealnej turbiny kocioł pracuje w sta- 
łych warunkach, czyli pochodna po czasie energii zgromadzonej w kotle jest zerowa. 


d 
g Er) (IV.18) 


Bilans pochodnej po czasie energii kotła uwzględnia przepływ 7 moli gazu na sekundę 
w analogiczny sposób jak przy przepływie gazu przez turbinę, ale zamiast wykonanej pracy 
uwzględnia moc cieplną, — z jaką ogrzewany jest gaz w kotle, oraz pomija wszelkie 
straty ciepła: 


d ; : 
0 z dt (Eociot ) =n (Win + Pinp Vinp ) = (U w Pout Vout » F O sana (IV. 19) 


Korzystajac z wprowadzonego pojecia entalpii, dla dowolnego gazu moc kotla mozna 
obliczyé: 


Ò aj m ń ` Ah (IV.20) 


Natomiast moc kotła podgrzewajacego gaz doskonały wynosi: 


jem u ń i CAT (IV.21) 


Oczywiście moc chłodzenia dla idealnej chłodnicy opisana jest identyczną zależnością. 
Idealny kocioł będzie się różnił od chłodnicy jedynie tym, że moc grzania Oi kotta jest 
wielkością dodatnią, a chłodnicy — ujemną. 

Ze schematu idealnego kotła widać, że ciśnienia na wejściu i wyjściu muszą być takie 
same p, =Po Oznacza to, że wykonanie bilansu pochodnej entropii dla idealnego kotła, co 
byłoby bardzo trudne, o ile wykonalne, na szczęście nie jest konieczne, gdyż nie prowadzi do 
uzyskania żadnej nowej informacji. 
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IV.3. Cykl Braytona: 
kociol, turbina, chlodnica i sprezarka 


Zostanie rozważony zamknięty cykl Braytona. Cykl ten jak dotąd nie znalazł zastoso- 
wania w energetyce, ale możliwe, że kiedyś zostanie zrealizowany w elektrowni jądrowej 
wykorzystującej wysokotemperaturowy reaktor jądrowy chłodzony helem — high temperatu- 
re reactor, w skrócie HTR. 


sprężarka 


Rysunek IV.3. Elektrownia gazowa z zamkniętym cyklem Braytona, 
w której wysokotemperaturowy reaktor jądrowy typu HTR chłodzony helem pełni funkcję kotła 


W przedstawionym na rysunku IV.3 układzie reaktor HTR ma pełnić funkcję kotła, 
a pozostałe jego elementy to turbina, chłodnica i sprężarka. Ważne też, aby tak dobrać 
parametry pracy turbiny, sprężarki i generatora, by były one połączone jednym wałem 
i moc, z jaką pracuje turbina, służyła do odpowiedniego napędzania sprężarki i genera- 
tora prądu. 


kocioł 
(reaktor HTR) 


or 


cisnienie p 


A chłodnica D 


objętość molowa v 


Rysunek [V.4. Zamknięty cykl Braytona w przestrzeni p-v 
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Na rysunku IV.4 można prześledzić w przestrzeni p-v parametry jednego mola gazu 
w cyklu, w jakim obiega on układ pokazany na rysunku IV.3. Kolejne procesy to: 

AB — adiabatyczne sprężanie — sprężarka, 

BC — izobaryczny ogrzewanie — kocioł, 

CD = adiabatyczne rozprężanie — turbina, 

DA = izobaryczne chłodzenie — chłodnica. 


Uwaga 1: W przypadku maszyn przepływowych wskazanie, że na osi X pokazana 
jest objętość gazu zamiast objętość molowa gazu przepływającego przez układ, 
może prowadzić do nieporozumień i pomyłek. 


Uwaga 2: Typowa elektrownia gazowa pracuje w otwartym cyklu Braytona, gdzie 
rolę chłodnicy pełni system pobierania z atmosfery zimnego powietrza, które wraz 
z gazem ziemnym trafia do komory spalania, a po spaleniu do atmosfery odprowa- 
dzane są ciepłe spaliny. Typowa turbina jest zintegrowana ze sprężarką i kotłem. 
W szczególności mówi się o komorze spalania w turbinie. 


Przedstawiony zarys opisu maszyn przepływowych pokazuje, że entalpia jest bardzo 
użytecznym pojęciem. Moc idealnej, adiabatycznej turbiny i sprężarki oraz moc idealnego 
kotła i chłodnicy są równe iloczynowi zmiany entalpii molowej przepływającego gazu i ilo- 
ści moli gazu przepływającego przez każdą z tych maszyn. 

Co więcej, aby obliczyć zmianę entalpii gazu, nie potrzeba znać szczegółów technicz- 
nych sposobu przechodzenia gazu ze stanu początkowego do końcowego. Wystarczy znać 
parametry początkowe i końcowe gazu. Dlatego mówi się, że entalpia jest potencjałem ter- 
modynamicznym, podobnie jak poznana energia wewnętrzna i entropia. Warto też podkre- 
ślić, że dla gazu doskonałego obliczenie zmiany entalpii molowej (entalpii jednego mola 
gazu) jest wyjątkowo proste: Ah = C AT. Analogicznie entalpia n moli gazu doskonałego 
wynosi: AH=n: Ah=n: CAT. 

Pewna trudność polega na tym, że w zależności od funkcji maszyny przepływowej 
zmiana entalpii opisuje albo pracę mechaniczną jak w turbinie, albo dostarczone ciepło jak 
w kotle. Bez wątpienia taka sytuacja jest trudna do zaakceptowania z intuicyjnego punktu 
widzenia i rodzi pytania o interpretację pojęcia entalpii, o wyjaśnienie, czym jest entalpia. 
Wydaje się, że najlepszym sposobem rozwiązania tej trudności jest zyskanie umiejętności 
rozwiązywania zadań, w których pojęcie entalpii jest użyteczne. Dzięki temu można nabrać 
pewności, że entalpia w odwracalnych maszynach przepływowych jest czymś w rodzaju 
waluty, wskazującej, ile należy zapłacić za dostarczone ciepło albo za wykonaną pracę. 


IV.4. Entalpia w nieodwracalnym procesie 
Joule’a—Thomsona (JT) 


Rysunek IV.5 ilustruje przebieg doświadczenia Joule’a-Thomsona (JT), w którym gaz 
przepycha się przez porowaty korek. Rury doprowadzające i odprowadzające gaz są ide- 
alnie izolowane termicznie, a porowaty korek jest nieważki (bardzo lekki), czyli nie może 
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pochłaniać ciepła, jego pojemność cieplna jest niemal zerowa, pomijalnie mała, nie moze ani 
zwiększać, ani zmniejszać swojej temperatury. 


n [mol/s] 
Pinp Vinp Pout Vout 
Tinp Tout 


Zerowa wymiana ciepta 
pomiedzy otoczeniem u 
i układem rurka + 

porowaty korek 


Rysunek IV.5. Schemat doświadczenia Joule’a~Thomsona 


Wielkość porów w korku i ciśnienia gazów są tak dobrane, aby przepływ gazu był 
kwazistacjonarny, by po obu stronach korka gaz był w równowadze termodynamicznej, 
czyli w szczególności miał dobrze określoną temperaturę i ciśnienie. Taki przepływ gazu 
jest dość powolny, równowagowy, ale w oczywisty sposób jest to proces nieodwracal- 
ny, gdyż przepływ gazu jest możliwy jedynie w kierunku od ciśnienia większego Pi do 
ciśnienia mniejszego p,,. Innymi słowy, jeśli nakręcić film z przebiegu doświadczenia JT 
i wyświetlić go od końca, to każdy zauważy, że w rzeczywistości jest niemożliwe, aby 
porowaty korek mógł spowodować przepływ gazu od ciśnienia mniejszego do większego. 

Energia korka jest pomijalnie mała, a energia gazu w korku Blue musi być stabilna, 


czyli pochodna po czasie energii zgromadzonej w korku musi być zerowa: 


d 
ae (E se torek ) = 0 (IV. 22) 


Zatem bilans pochodnej po czasie z energii gazu w korku jest następujący: 


d . 
0 = dt (E paz _ toret ) = (Cin + Pinp Vinp ) ~ (U + Pou Vout » (IV.23) 
Korzystając z wprowadzonego pojęcia entalpii, można zauważyć, ze w doświadczeniu JT 
dla dowolnego gazu zawsze zachodzi równość: 


out ig BRO (IV.24) 
Czyli w nieodwracalnym procesie JT entalpia przepływającego gazu się nie zmienia. 
Ten wynik bez wątpienia kłóci się z naszą intuicją, gdyż czujemy, że gaz o większym ciśnie- 
niu ma większą użyteczność niż gaz o ciśnieniu mniejszym. Jednak już wstępne zapozna- 
nie się z pojęciem entropii pokazuje, że o użyteczności gazu w maszynach przepływowych 
decyduje nie tylko entalpia, ale również jego entropia. Im entropia gazu jest mniejsza, tym 
użyteczność gazu jest większa. Wydaje się, że najlepszym sposobem opanowania trudnej 
umiejętności oceny użyteczności gazu jest nabranie wprawy w rozwiązywaniu zadań. 
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IV.5. Gaz doskonaly i gazy rzeczywiste 
w doświadczeniu JT 


Dla gazu doskonałego doświadczenie Joule’a—Thomsona jest zwyczajnie nudne, gdyż — 
jak pokazano — zmiana entalpii molowej gazu w doświadczeniu JT jest zerowa. Wiadomo 
również, że dla gazu doskonałego Ah = C AT, czyli jeśli Ah = 0, to i AT= 0. Coz może być 
ciekawego, gdy z obu stron korka temperatura gazu jest taka sama? 

Równocześnie widać, że pomiar AT w doświadczeniu JT dla gazów rzeczywistych 
może dostarczyć informacji o istotnej różnicy pomiędzy gazami rzeczywistymi i gazem 
doskonałym. Dla takich gazów jak azot i tlen doświadczenie Joule'a-Thomsona wyko- 
nane w temperaturze pokojowej pokazało, że po przeciśnięciu przez porowaty korek ich 
temperatura się obniża. Zauważyli to James Prescott Joule i William Thomson (Lord 
Kelvin) już w 1852 roku, ale nie przekuli tego wyniku na praktyczne zastosowania. Zro- 
bili to w Krakowie ponad dwadzieścia lat później, w 1883 roku, profesorowie Zygmunt 
Wróblewski i Karol Olszewski, skraplając azot i tlen. Zasada działania ich aparatury jest 
prosta. Gaz po przeciśnięciu przez porowaty korek zmniejsza swoją temperaturę i zostaje 
wykorzystany do schładzania gazu, który zmierza w kierunku porowatego korka i dopie- 
ro po wykonaniu tej misji jest wypuszczany do atmosfery. Dzięki temu kolejne porcje 
gazu zmierzające do korka mają coraz niższą temperaturę. Jeśli układ jest dobrze izolo- 
wany termicznie i do dyspozycji jest odpowiednio duży zbiornik ze sprężonym gazem, to 
widać, że prędzej czy później opisany mechanizm chłodzenia doprowadzi do tak niskiej 
temperatury, że gaz się skropli. Jeśli rura zostanie ustawiona pionowo, gaz będzie się 
przeciskał przez porowaty korek z góry na dół i po jakimś czasie zacznie z niego kapać 
skroplony gaz. 

Warto dodać, że w doświadczeniu JT przeprowadzanym w temperaturze pokojo- 
wej wodór i hel zwiększają swoją temperaturę. Dzieje się tak również w niższych tem- 
peraturach, co powoduje, że skroplenie szczególnie helu jest bardzo trudne i wymaga 
albo zastosowania (po wstępnym ochłodzeniu) odpowiedniego wysokiego ciśnienia, 
jak w skraplarkach Simona, albo też sprężenia i wstępnego ochłodzenia za pomocą 
ciekłego wodoru. 


IV.6. Paradoks 
w schładzaniu się gazu doskonałego 


W codziennych doświadczeniach spotykamy się z sytuacją, gdy czujemy, że zbiornik ze 
sprężonym gazem schładza się, kiedy gaz jest z niego wypuszczany. Tak dzieje się na przy- 
kład z pojemnikiem z dezodorantem w spreju. W wielu źródłach ten powszechnie odczuwany 
chłód pojemnika przywołuje się jako przykład ilustrujący omówione powyżej doświadczenia 
Joule'a-Thomsona. Podkreśla się, że obniżenie temperatury zbiornika jest możliwe jedynie 
dzięki temu, że wypuszczany jest z niego gaz rzeczywisty, a nie doskonały. 

Niestety taka interpretacja tego doświadczenia jest fałszywa, co ilustruje zadanie roz- 
wiązane w kolejnym rozdziale. 
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IV.7. Rozwiązanie paradoksu 
w schładzaniu się gazu doskonałego 


W celu uproszczenia zadania butla ze sprężonym gazem doskonałym zostanie unie- 
ruchomiona i umieszczona w nieskończenie dużej próżni, na przykład kosmicznej. Układ 
butli z gazem jest układem zamkniętym i izolowanym. Jeśli przez niewielki otwór gaz 
będzie wylatywał do próżni, to można użyć bilansu energii do opisu zachodzących pro- 
cesów. Jak dotąd wszystkie omawiane procesy były równowagowe, czyli wypadkowa siła 
wywierana na tłok była zerowa, a energia kinetyczna przepływających paczek z gazem 
była pomijana. Tu jednak należy uwzględnić niezerową siłę wypadkową, z jaką wypychana 
jest każda paczka z gazem i tym samym energię kinetyczną, z jaką paczka jako całość się 
porusza. Bilans energii przed wypuszczeniem i po wypuszczeniu ze zbiornika niewielkiej 
części mola gazu On jest następujący: 

nC T=(n-ón)C (T—67) +ón:CT+E 


u (IV.25) 
gdzie: 
nCT -energia wewnętrzna gazu w butli przed wypuszczeniem nie- 
wielkiej porcji gazu zawierającej ön moli, 
(n- ón)C (T—6T) — energia wewnętrzna gazu w butli po wypuszczeniu ön moli 
gazu, 
ön: CT — energia wewnętrzna paczki gazu, która opuszcza butlę, 
— energia kinetyczna, z jaką porusza się jako całość paczka 
wypuszczonego gazu. 


kin_6n 


Oczywiście paczka gazu wypuszczonego do próżni przestanie być paczką o ustalonej 
objętości, rozpłynie się w próżni, ale nadal centrum masy tej paczki będzie poruszać się ze 
stałą prędkością względem butli. Paczka ta wciąż będzie miała jako całość energię kinetycz- 
ną równą pracy, jaką wykonał gaz, wypychając ją z butli, czyli: 
= pov (IV.26) 


ZK 


gdzie SV jest objętością, jaką zajmowała w butli paczka licząca ön moli gazu tuż przed 
wypchnięciem jej z butli. Gaz w butli jest gazem doskonałym, czyli równanie stanu dla ön 
moli zajmujących objętość V jest też równaniem Clapeyrona: 


p: óV=6n:RT (IV.27) 
Zatem równanie (IV.25) bilansujące energię przybiera postać: 
nC T= (n— 0n)C (T-T) + õn: CT+6n: RT (IV.28) 


Po prostych przeksztatceniach i pominięciu wielkości (ön * 67) jako znacznie mniejszej 
niż wszystkie inne dostaje się: 


nC,:8T=ó8n-RT (IV.29) 
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Układając równanie różnicowe, przyjęto, że zmiana temperatury ôT i liczba moli wyrzu- 
canego gazu ön są wielkościami dodatnimi. Jednak zarówno temperatura gazu w butli, jak 
i liczba moli gazu w butli zmniejszają się. Zatem przechodząc od równania różnicowego do 
różniczkowego, należy dokonać konwersji znaków: 


970 -dT, ön —> -dn (IV.30) 


Otrzymujemy w ten sposób równanie różniczkowe: 


nC dT = RTdn (IV.31) 


Równanie to łatwo rozwiązać, wykorzystując stosowaną już w tej publikacji metodę 
rozdzielenia zmiennych, czyli zapisując je w takiej postaci, aby każdy jego człon zależał 
tylko od jednej zmiennej: 


dr R dn 


TO n (IV.32) 


Widać, ze obie strony można scałkować od stanu początkowego do końcowego. Po 
lewej stronie stan początkowy i końcowy opisują temperatury 7, i 7,. Po prawej stronie rów- 
nania zaś — początkowa i końcowa liczba moli gazu w butli, odpowiednio n, i n,. Całkowe 
równanie ma zatem postać: 


| = Í = IV.33 
Ty T C, ng n ( ) 
Jego rozwiązaniem jest: 
de R ny 
h| = |=—- Inj) = IV.34 
T ) C, l m | ? 
co po prostym przekształceniu daje: 
R 
n, \& 
T=] = (IV.35) 
no 


Dla jednoatomowego gazu doskonałego, C, = (3/2)R, gdy w zbiorniku pozostanie połowa 
gazu, n, — n,/2, jego temperatura końcowa wyniesie 7, ~ 0,63 : T, 

Podsumowując, pokazano, że gdy gaz doskonały jest wypuszczany do próżni, to schła- 
dza się butla, w której gaz ten był zgromadzony. Dla przypomnienia: w przypadku rozprę- 
żenia gazu do drugiego zbiornika” energia każdej paczki wypuszczonego gazu pozostawała 
w układzie dwóch izolowanych zbiorników i temperatura układu pozostawała stała. 

Zadanie podobne do przedstawionego bardzo często pojawia się na egzaminach wstęp- 
nych na studia doktoranckie na najlepszych uczelniach technicznych w USA, a także w Polsce. 


Patrz: rozdział 4.1: Procesy odwracalne i nieodwracalne, w: Fizyka dla szkół wyższych. Tom 2, 
dz. cyt. 
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V. Termodynamika ogniwa paliwowego 


V.1. Spalanie wodoru w tlenie 
w przepływowej komorze spalania 


Na początku zostanie przeanalizowany proces całkowitego spalania wodoru w tlenie 
w przepływowej komorze spalania pokazanej na rysunku V.1. System jest tak skonstruowa- 
ny, aby ciśnienie w całym układzie było stałe p = const, stała była temperatura pracy komory 
spalania T = const i była równa temperaturze wodoru, tlenu oraz powstającej wody. Stałość 
temperatury zapewnia odpowiedni system chłodzenia komory spalania. 


T = const, p = const 


Rysunek V.1. Przepływowa komora spalania wodory w tlenie 
w warunkach stałej temperatury i stałego ciśnienia 


Każdy mol wodoru i tlenu, wpadając do komory spalania, ma ciśnienie p,, i po, obję- 
tość molową (czyli objętość jednego mola) v,, iv, , oraz temperatury 7,, i T, a opuszcza- 
jaca komorę woda (gaz, para wodna) ma parametry p, o Vy os Tmo: Spełnione są też warunki 
Pr, Po Pia “Peons OTAZ Ta, = To, = Tiyo lt 

Należy podkreślić, że podobnie jak w przypadku turbiny i kotła pomijana jest energia 
kinetyczna paczek gazu przepływających przez komorę spalania. Innymi słowy zakłada się 
zerową prędkość unoszenia, z jaką gazy przepływają przez komorę spalania. 

Z ogólnej wiedzy i ze schematu pokazanego na rysunku V.1 wynika, że idealna komo- 
ra spalania, pracując w stabilnych warunkach, musi byé odpowiednio chłodzona, czyli 
moc O ciepła dostarczanego do komory ma wartość ujemną. Aby komora pracowała 


grzenie 


H0? 
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stabilnie, energia zgromadzona w komorze spalania E 


Komora MUSİ być stała w czasie, czyli jej 
A . 3 omora 
pochodna po czasie musi być zerowa: 


d 
F (E morn) = 9 (V.1) 


Całkowite spalenie wodoru w tlenie ma szanse zajść (i zachodzi), gdy do komory wpada 
dwa razy mniej moli tlenu niż wodoru i gdy spełniona jest stechiometryczna proporcja liczby 
moli gazów: 


1 
1-H, +>-0, >1-H,O (V.2) 
Kazdy mol gazu przepływający przez komorę zmienia jej energię o wartość równą ental- 


pii molowej A = u+ pv i widać, że tlen oraz wodór zwiększają energię komory, a opuszcza- 
jąca ją woda — zmniejsza. Zatem warunek stacjonarności pracy komory ma postać: 


d ; 1 . 
0= dt (E omora ) EG (h + 2 ho, E No ) + Ooms (V.3) 
gdzie hi» h öp ho to entalpie molowe wodoru, tlenu i wody. Z równania (V.3) wprost dosta- 
je sie: 
‘ : 1 
Oi =n hr o ~ hy, + 7% (V.4) 
Og. = ń í (h a hi») (V.5) 


gdzie w równaniu (V.5) wprowadzono oznaczenie na wejściową i wyjściową entalpię molo- 
wą: hz =h, + (U/2)h,, oraz h = Aao: Zostanie teraz wprowadzona wielkość tymczasowo 
nazwana zmianą entalpii molowej Ah = ha hing: Widać, że inaczej niż w przypadku turbiny 
i kotła zmiana entalpii uwzględnia efekt zniknięcia wodoru i tlenu oraz stworzenia wody. 
Potrzebna jest lepsza nazwa i powszechnie wielkość ta nazywana jest entalpią tworzenia, 
tu entalpią tworzenia mola wody. Standardowa wartość entalpii tworzenia wody, czyli przy 
stałej temperaturze T = 298 K i stałym ciśnieniu p = 1 bar = 10° Pa, dla wody w postaci 
gazowej wynosi: 


out 


= —241,60 kJ/mol, 


0 = 
298 (g H20) 


a w postaci ciekłej, czyli po skropleniu: 


AH? = —285,84 kJ/mol. 


298 (c H>0) 


Wielkości liczbowe entalpii tworzenia mola wody są podane przy założeniu, że standar- 
dowe entalpie tworzenia wodoru i tlenu w postaci gazowej są dokładnie zerowe: 


AM yg (gH) 0,00 kJ/mol, 
AH? go, = 0,00 kJ/mol. 
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Podsumowując, moc cieplna wydzielana przy całkowitym spalaniu wodoru w tlenie, 
gdy produkowana woda jest w postaci ciekłej, wynosi: 


O dagen: = -Ò „ania (V. 6) 
O pasie =h: 285,84 kJ/mol 


Z przeprowadzonej analizy wynika, że w reakcjach egzotermicznych entalpia tworze- 
nia jest ujemna. Intuicja podpowiada też, że opisany proces całkowitego spalania wodoru 
w tlenie jest procesem nieodwracalnym. Film nakręcony podczas takiego procesu i pusz- 
czony od końca pokaże — wydaje się niemożliwą do realizacji — sytuację, gdy para wodna 
podgrzewana w komorze rozkłada sie na wodór i tlen, i w dodatku te dwa gazy zostają 
rozdzielone. 


V.2. Odwracalna reakcja wodoru z tlenem 
w przepływowej komorze reakcji 


Dysponując analizą zmian entropii, można zaprojektować maszynę przepływową 
spełniającą termodynamiczny warunek odwracalności. Maszyna taka, zasilana parą wod- 
ną i energią elektryczną, będzie rozkładała wodę na wodór i tlen. W odwrotnej sytuacji, 
czyli gdy będzie zasilana wodorem i tlenem, gazy te, łącząc się w cząsteczki wody, będą 
generowały energię elektryczną. Pierwszy typ maszyn to elektrolizer, który pobiera energię 
elektryczną i ciepło, a drugi to ogniwo paliwowe generujące energię elektryczną, ale wyma- 
gające też chłodzenia. Należy podkreślić, że techniczne szczegóły każdej z tych maszyn 
przepływowych muszą być różne i na rysunku V.2 są intencjonalnie całkowicie pominięte, 
aby podkreślić, że przeprowadzona poniżej analiza określi jedynie termodynamiczne warun- 
ki odwracalności. 


a) = const, = const 5 T = const, p = const 


i LE (pe 


I A Hz | 
Ran =i H20 

= i—> | — 
=P 


—— 


odbiór <— | zasilanie 
energii I energią 
lekt lekt. 
elektrycznej 02 elektryczną 
ffit o ae w] tt o Ogszante | w] 
dodatnie 


Rysunek V.2. Idealne ogniwo paliwowe (a) i idealny elektrolizer (b) 
jako dwie idealne, wzajemnie odwracalne maszyny przepływowe 
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V.3. Ogniwo paliwowe 


Dalsze rozważania zostaną przeprowadzone dla idealnego, odwracalnego ogniwa pali- 
wowego. Będzie można je wprost wykorzystać do analizy pracy elektrolizera, gdyż obie te 
maszyny są wzajemnie do siebie odwracalne. Warunek stacjonarności pracy ogniwa paliwo- 
wego można zapisać jako: 


d f 1 . 
0 = dt (E ogniwa ) =N; (ba + 2 ho, s hyo + O rania z ra (V.7) 


gdzie P.„„, to moc elektryczna wyprowadzona na zewnątrz ogniwa. Znak plus przed Ò „sania 
i znak minus przed P.„„,„ odzwierciedlają przyjętą konwencję pierwszej zasady termodyna- 
miki, AU = Q — W, gdzie Q oznacza ciepło dostarczone do układu (tu: ogniwa paliwowego), 
a W pracę wykonaną przez układ na zewnątrz tego układu, tu: wyprowadzoną na zewnątrz 
moc elektryczną P in 

Warunek odwracalności oznacza, że dostarczona z zewnątrz moc cieplna musi być 


związana z dostarczaną z zewnątrz entropią przez relację: 


Oe. 
AS kawie = (V.8) 
z zewnatrz Toorat 
gdzie T „Jest stałą temperaturą pracy ogniwa paliwowego. Po drugie, entropia ogniwa pali- 
wowego także musi być stała, czyli spełniać warunek: 
d i 1 . 
0 = dt Gaj =ne SH, at 5% z Si,0 de AS postarczone (V.9) 


z zewnątrz 


gdzie s to entropie molowe wodoru, tlenu i wody. 


H? Soy Sino 


Korzystając z równań (V.7)-(V.9), dostaje się moc elektryczną ogniwa: 


i 1 
Poore = H+ (o m dami i Sq,o) u G m Ii ` SH, ) + 2 (ho, m dus ` So, ))} 


Widać, że warto zdefiniować nową funkcję stanu: 


def 


g=u=h-Ts=u+ pv-Ts (V.10) 


Czasem jest ona oznaczana małą litera g i nazywana molową entalpia swobodną, molo- 
wym potencjałem Gibbsa, ale najczęściej jest nazywana potencjałem chemicznym i oznacza- 
na grecką litera u. Po pomnożeniu przez liczbę moli dostaje się funkcję stanu oznaczaną dużą 
literą G i nazywaną entalpią swobodną lub potencjałem Gibbsa: 


def 
G=H-TS=U+pV-TS (V.11) 
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Dygresja: Warto dodać, że wprowadzając termodynamike na podstawie spójnego 
formalizmu matematycznego, najpierw definiuje się potencjał Gibbsa G — równanie 
(V.11), a następnie potencjał chemiczny każdego ze składników mieszaniny jako 
pochodną cząstkową potencjału Gibbsa po liczbie moli tego składnika n, przy stałym 
ciśnieniu, stałej temperaturze i stałej liczbie moli wszystkich pozostałych składników 


def 
u, =(0G /0h, ), „„,,. Dalej dowodzi się, że G = Ż tlle Koniec dygresji. 


Po wprowadzeniu pojęcia molowej entalpii swobodnej moc elektryczną można teraz 
zapisać w bardziej przejrzysty sposób: 


: 1 
Piy “A (o [eu + 7 80: ) (V.12) 


A seket (Sai m 5, 


A następnie, podobnie jak w przypadku entalpii, można wprowadzić pojęcie molowej 
swobodnej entalpii tworzenia AG oraz standardowej molowej entalpii swobodnej tworze- 
nia. Wartość standardowej swobodnej entalpii tworzenia, czyli przy stałej temperaturze 
T = 298 K i stałym ciśnieniu p = 1 bar = 10° Pa, dla wody w postaci gazowej wynosi: 


AG? = —228,57 kmol, 


298 (g H30) 


a w postaci ciekłej, czyli po skropleniu: 


AG? =—237,13 kJ/mol. 


298 (c H20) 
Wielkości te pochodzą z pomiarów i dają wprost oszacowanie mocy elektrycznej ideal- 
nego ogniwa paliwowego, gdy z ogniwa odbierana jest woda w postaci ciekłej. 


P an 711 237,13 kJ/mol (V.13) 
Pamiętając, ze moc cieplna wydzielana przy całkowitym spalaniu wodoru w tlenie (V.6) 
wynosi ELMENT = n 285,84 kJ/mol, można określić sprawność termodynamiczną ogniwa 
paliwowego: 
Pigs n- 237,13 
= = — = 83% 
ý wdówka n- 285, 84 : ce Ę 


Widać, że sprawność idealnie odwracalnego ogniwa paliwowego pracującego w stałej 
temperaturze T = 298 K i pod stałym ciśnieniem p = 1 bar = 10° Pa jest bardzo wysoka, co bez 
wątpienia jest inspiracją do konstruowania coraz lepszych ogniw paliwowych. Należy dodać, 
że przy tak niskiej temperaturze reakcje zachodzą dość powoli i aby zwiększyć moc, konstru- 
owane są wysokotemperaturowe ogniwa paliwowe pracujące przy temperaturach kilkuset 
stopni Celsjusza. Wówczas jednak ich sprawność maleje. Co więcej, całkowicie pominięte 
tu rozwiązania konstrukcyjne, w tym dobór odpowiednich materiałów, z których wytwarzane 
są elektrody oraz elektrolit, nadal jest wyzwaniem i dlatego ogniwa paliwowe nie są jeszcze 
powszechnie wykorzystywane. 
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VI. Warunek spontanicznego 
przebiegu reakcji chemicznej 


Na przykładzie ogniwa paliwowego pokazano, że spontaniczność reakcji chemicznej 
w warunkach stałej temperatury i stałego ciśnienia związana jest z możliwością generowania 
energii elektrycznej, co pokazuje ujemna wartość zmiany entalpii swobodnej (potencjału 
Gibbsa) AG < 0. Odwrotny przebieg reakcji jest możliwy, ale wymaga dostarczenia ener- 
gii elektrycznej do elektrolizera, co pokazuje AG > 0. Oczywiście chemia daleko wykracza 
poza problem ogniw paliwowych i elektrolizerów, a wprowadzona entalpia swobodna, którą 
można określić dla każdej reakcji chemicznej, dostarcza informacji o tym, czy dana reakcja 
będzie przebiegać spontanicznie w warunkach ustalonej temperatury i stałego ciśnienia: 


Jeśli zmiana entalpii swobodnej jest ujemna, AG < 0, to przy zapewnieniu sta- 
łej temperatury i stałego ciśnienia reakcja przebiega spontanicznie. 


VII. Synteza amoniaku 


Już w XIX wieku odkryto, że nawozy azotowe znacznie zwiększają plony, ale nie potra- 
fiono ich syntetyzować. Niemal jedynym ich źródłem była saletra chilijska, będąca również 
surowcem do produkcji materiałów wybuchowych, a utrudniony dostęp do niej ograniczał 
zbrojenia. Niektórzy stawiają tezę, że pierwsza wojna światowa mogła wybuchnąć jedynie 
dlatego, że Niemcy, a za nimi inni, opanowali technologię syntezy amoniaku. Szacuje się, 
że we współczesnym Świecie dzięki nawozom sztucznym zaspokojony jest głód miliardów 
ludzi. Szacuje się też, że materiały wybuchowe powstałe na bazie amoniaku i technologii jego 
syntezy spowodowały śmierć nawet miliarda ludzi. Technologia syntezy amoniaku pokaza- 
ła też drogę do syntezy innych produktów, w tym gazów bojowych, takich jak cyklon B, 
którym Niemcy wymordowali miliony Żydów. Bez wątpienia ludzkość musi pracować nad 
odpowiedzią na pytania o znaczenie etyki w nauce. Niemniej ważne jest też, aby wiedzieć, 
jak zsyntetyzować amoniak. 

Reakcja spalania wodoru w tlenie w stałej temperaturze T = 298 K i pod stałym ciśnie- 
niem p = 1 bar przebiega spontanicznie. Co więcej, spalenie wodoru jest praktycznie całko- 
wite. Analogicznie jak w przypadku spalania wodoru w tlenie szkolna chemia podpowiada, 
że aby zsyntetyzować mol amoniaku, reakcja musi przebiegać według następującej proporcji 
liczby moli: 

3 1 
> + > N, > NH, (VII.1) 

O kierunku przebiegu spontanicznej reakcji syntezy lub rozpadu amoniaku decy- 
duje znak swobodnej entalpii tworzenia amoniaku i wiadomo, ze w standardowym stanie 
(T= 298 K, p= | bar) dla syntezy mola amoniaku AGH, gut) ~~ 16,45 kJ/mol. Można zatem 
oczekiwać, że synteza zajdzie spontanicznie, ale niestety tak nie jest. Jednym z powodów 
jest stosunkowo mała ujemność swobodnej molowej entalpii tworzenia. Gdy reakcja osią- 
gnie stan równowagi, czyli gdy liczba reakcji syntezy będzie równa liczbie reakcji rozpadu 
amoniaku, to produktem nie będzie czysty amoniak, ale mieszanina trzech gazów. Co więcej, 
w tej mieszaninie praktycznie nie będzie amoniaku. Nasuwa się zatem pytanie, w jakich 
warunkach, czyli w jakiej temperaturze i pod jakim ciśnieniem, należy przeprowadzić reak- 
cję, aby produkować amoniak. 

Można wykonać serie doświadczeń dla różnych stałych temperatur, a dla kolejnych sta- 
łych temperatur dla różnych ciśnień. Po ich wykonaniu uzyska się wiedzę, dla jakiej tempe- 
ratury T 1 dla jakiego ciśnienia p, produkcja amoniaku będzie najwydajniejsza. Warto 


const const 


jednak najpierw wykonać poniższą analizę. 
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Podana powyżej wartość AG", wi, 71645 kmol odnosi się do teoretycznej sytua- 
cji całkowitego przereagowania wodoru i azotu w amoniak, co ilustruje rysunek VII.1. Tuż 
przed reakcją gazy A i B (odpowiednio wodór i azot) znajdują się w dwóch stykających się 
idealnych cylindrach z ruchomymi tłokami. W każdym ze zbiorników jest zapewniona stała 
temperatura odniesienia T, „, „= 298 K i stale ciśnienie odniesienia p „„, „= 1 bar. Po usunię- 
ciu przegrody dzielącej zbiorniki A i B oraz po całkowitym przereagowaniu gazów w jednym 
wspólnym zbiorniku będzie znajdował się jedynie gaz C (amoniak) i zostanie zapewniona 
jego stała temperatura T = 298 K i stale ciśnienie p = | bar. 


const 0 const 0 


Przed reakcja 


> > 
Pow Pconst_0 Pconst_0 = 1 bar Few 
Tconst_0 Teonst_0 = 298 K 
> > 
Frew Pconst_o = 1 bar | fay 
Gr | ————— 


Teonst_0 = 298 K 


Rysunek VII.1. Reakcja powstania produktu C po usunięciu przegrody dzielącej zbiorniki Ai B 
i całkowitym przereagowaniu substratów A i B w warunkach stałej temperatury i stałego ciśnienia 


Parametry charakteryzujące przebieg tak wyidealizowanej teoretycznej reakcji można 
zapisać w tabeli VII.1. 


Tabela VIL1 


Parametry całkowitego przereagowania wodoru z azotem w amoniak 


Gaz A Gaz B Gaz C AG 
(wodór) (azot) (amoniak) 
i : 3 1 
Stechiometria —H, =N, NH R 
2 2 3 
Liczba moli 3 1 0 
przed reakcja 2 2 kJ 
AG roan, = —16,45 aa 

Liczba moli mo. 

5 0 0 1 
po reakcji 


Realne reakcje chemiczne prawie nigdy jednak tak nie zachodzą i w stanie równowagi 
chemicznej zawsze będzie występowała mieszanina trzech gazów. Dla syntezy dwóch sub- 
stratów (A, B) w jeden produkt (C) w warunkach stałego ciśnienia p „„„ i stałej temperatury 


const 


T parametry charakteryzujące stan równowagi zapisano w tabeli VIL.2. 


const 
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Tabela VII.2 


Parametry stanu równowagi w reakcji syntezy produktu C z substratów A i B 


Substrat A Substrat B Produkt C Entalpia 
swobodna 
Stechiometria vA vB Ve C é 
Liczba moli przed reakcją Vy Vp 0 — 
Liczba moli a b c 
"S| Ciśnienia parcjalne p p P A 
y paroj 4 2 E Minimalna 
ô 25 
a b wartosé G 
z Wz = = ae 
2 a+b+c a+b+c a+b+c rowna 
s | Względna liczba moli C omona 
rn Xi = P A X; = P B X = P C 
Peonst Pronst P; onst 
Chciejstwo X, =0 X =0 Xo =l chciejstwo 


Dojście do stanu równowagi jest procesem nieodwracalnym, a gdy zostanie on osią- 
gnięty, to entalpia swobodna będzie najmniejsza. Niemożliwy jest dalszy spontaniczny 
postęp reakcji, niemożliwe jest osiągnięcie stanu chciejstwa, całkowitego przereagowania 
substratów A i B w produkt C. W stanie równowagi liczby moli każdego z gazów wyniosą 
odpowiednio a, b, c, względne liczby moli X, X,, X,. Suma ciśnień parcjalnych gazów A, B 
i C będzie równa stałemu ciśnieniu, przy którym zachodzi reakcja p, „=p, + p, p,. Stan 
równowagi chemicznej ilustruje rysunek VII.2. 


Stan równowagi mieszany gazów A, B, C 


> 
Fzew 


Pconst = Pat DB Pc 
Xa Xg, Xc I 


Tconst 


Rysunek VII.2. Reakcja powstania równowagowej mieszaniny A, B i C 
w warunkach stałej temperatury i stałego ciśnienia 


Entalpia swobodna jest wielkością ekstensywną, stanowi sumę entalpii swobodnych 
wszystkich trzech gazów. W stanie równowagi wartość absolutna, której nigdy nie poznamy, 
ale nie będzie nam to potrzebne, będzie wynosiła: 
+b: utc: We (VII.2) 


; 20S 
równowagi u A 


gdzie u, U, Ho to molowe entalpie swobodne każdego z gazów, które dla wygody war- 
to tu nazywać potencjałami chemicznymi gazów A, B, C. Każdy stan, w tym stan równo- 
wagi, uwzględnia fakt, że każdy z trzech gazów zajmuje całą dostępną objętość, że suma 
ciśnień parcjalnych tych gazów równa jest stałemu ciśnieniu, przy którym zachodzi reakcja: 
P=Pyt Pst Pe 
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W stanie nieznacznie oddalonym od stanu równowagi liczba moli poszczególnych 
gazów niewiele się zmieni zgodnie ze stechiometrią reakcji i proporcjonalnie do dowolnie 
małej, ale niezerowej liczby e. Liczba moli gazu C wzrośnie o ev_, a liczba moli gazów Ai B 
zmniejszy się odpowiednio o ev, i ev, Wartość entalpii swobodnej mieszaniny gazów będzie 
wówczas wynosiła: 

GO =(a—e'v)':u,t(b-e'v,)upt(cte'v.)'H, (VII.3) 


równowagi 


Skoro e jest dowolnie małą, ale niezerową liczbą, to stan AE również będzie sta- 
nem równowagi i różnica wartości entalpii swobodnych dla tych dwóch stanów będzie zero- 
wa, G® —G, „= 0. Po prostych rachunkach i podzieleniu równania stronami przez 

równowagi rownowagi 


e dostaje się warunek równowagi chemicznej: 


Volo — WU, + Vahl) = 0 (VII.4) 


W tym momencie dla każdej stałej temperatury T „, wprowadza się model, w którym 
zmiany potencjałów chemicznych u każdego z gazów będą obliczane tak, jakby były one 
gazami doskonałymi. Zasadność tego modelu opiera się na obserwacji, że przy mieszaniu 
gazów doskonałych entropia i entalpia swobodna bardzo silnie zależą od proporcji miesza- 
nych gazów. Dlatego można oczekiwać, że efekt mieszania będzie tak duży, że fakt, iż gazy 
A, B, C nie są gazami doskonałymi, a rzeczywistymi, niewiele zmieni. 

Zatem w modelu będzie badana zmiana potencjału chemicznego gazu (swobodnej 
entalpii molowej) dla ustalonej temperatury T „, w zależności od ciśnienia gazu. Dla gazu A: 


cons: 


AW ap, = ACA —T, 


A(Pa) const SA PA 


) = Aa.) = T ons ASA) (VII.5) 


Ale wiadomo, ze dla gazu doskonałego Ah = 0, gdy temperatura jest stała. Wiadomo 
również, że zmiana entropii molowej (patrz wzór (II. 14)) wyniesie: 


ść 
AS ap) =C, tof = | R: uf? 4 | =-R. (2) (VIL6) 


gdzie p, = 1 bar jest ciśnieniem odniesienia. 


Po podstawieniu zależności (VII.6) do równania (VII.5): 


Ab ap 1) ej O i m[ 24) (VII.7) 


Wielkość A ,(,,, można też zapisać jako różnicę: 
Akce) = Hace 7 HA (VIL8) 


gdzie u, to potencjał chemiczny dla stałej temperatury T „, i stałego ciśnienia odniesienia 


const 
p,= 1 bar. 
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Teraz można zapisać bezwzględne wartości potencjału chemicznego gazów A, BiC dla 
stałej temperatury 7, ,, w funkcji zmiany ciśnienia tych gazów jako: 


const 


Hp) = My ERI st g In 2) (VII.9) 
Po 


p 
Heeps) = W ERT cup ‘In z) 


P 
HC) = ue ERI oi ‘In a 


gdzie p , p,, p, to ciśnienia parcjalne gazów A, Bi C wstanie równowagi, gdy stałe ciśnienie, 
przy którym zachodzi reakcja, wynosi: p = p, + p, + p,. Pamiętając, że w stanie równowagi 
zachodzi równość (VII.4) i po podstawieniu do tego równania trzech równań (VII.9), dostaje 
się dla ustalonej temperatury T 
(2) 
E gine SE ONE 


const" 
0 = [vene RZY +v) | const va VB 
Pa Po 


A G®Teonst = RATA = s, + ven] (VII.11) 


3 | 
Po 
K, = Vy VB 
(2) (2) 
Po b, 


co pozwala w skrócie zapisać warunek równowagi jako: 


(VII.10) 


Dalej wprowadza się oznaczenia: 


0 = AG” const + RT ln(K,) (VII.12) 


const 


Wielkość AG jest określona dla ciśnienia referencyjnego p, = | bar i nie zależy od 
ciśnienia, a to oznacza, że wielkość K, również nie może zależeć od ciśnienia, które jest sumą 
ciśnień parcjalnych p =p, + p, + Po Znajomość ciśnień parcjalnych umożliwia wyznaczenie 
względnego udziału każdego z gazów w mieszaninie: 


X =, X, =, X, = (VIL.13) 
P P P 
Podstawiając: 
Pa _ Pa P, Po _ Pa P, a a 
Po P Po Po P Po Po P Po 
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K można zapisać jako: 
p 


ve-(vatva) 8 
xv 
K,= (2) (| (VIL.15) 
Po XX; 


Warto też wprowadzić wielkość K „zależną jedynie od stechiometrii reakcji i od względ- 
nej liczby moli gazów w stanie równowagi: 


X 
Ky = a (VII.16) 
x = Xv 


Ma ona tym większą wartość, im w mieszaninie w stanie równowagi jest więcej gazu C. 
Teraz równanie (VII.15) można zapisać jako: 


c-(vatva) 
K,= (2) -Ky (VII.17) 
Po 


Skoro K nie zalezy od cisnienia, to tym wieksze bedzie stezenie gazu C, im wielkosé 
t= (pipe °4*’» będzie mniejsza. W przypadku syntezy amoniaku v, = 3/2, v, = 1/2,v.= 1 
i wielkość t zależy od ciśnienia w następujący sposób: t = (p/p,)', czyli im większe będzie 
ciśnienie, tym t będzie mniejsze i tym większe będzie stężenie amoniaku w stanie rownowa- 
gi w mieszaninie wodoru, azotu oraz amoniaku. 

Po raz pierwszy wysokociśnieniową syntezę amoniaku przeprowadził Fritz Haber — 
w 1908 roku, uzyskując kilkuprocentowe stężenie tego gazu. 


VIII. Zakonczenie 


Analiza pokazana w niniejszej publikacji jest bazą wykładów z termodynamiki tech- 
nicznej i termodynamiki chemicznej, gdzie najczęściej milcząco zakłada się, że wszystkie 
potencjały termodynamiczne (funkcje stanu), w tym entropia, entalpia i entalpia swobodna, 
są już dobrze znane. Autor ma nadzieję, że informacje zawarte w książce przyczynią się do 
lepszego poznania i zrozumienia tych pojęć już na etapie nauki podstaw termodynamiki. 
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